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Entanglement and Ground States of Gapped Hamiltonians 

Abstract 

This thesis weaves together three separate results whose common thread is the 
study of quantum entanglement, a physical resource which, much like energy, may be 
exploited to perform tasks with the potential for extraordinary applications in areas 
of security, computing and simulation of classically intractable quantum phenomena. 

We begin by considering entanglement properties of an important class of quan- 
tum states, introduced by Fannes, Nachtergaele and Werner, known as Finitely Cor- 
related States (FCS). We derive bounds for the entanglement of a spin with an (adja- 
cent and non-adjacent) interval of spins in an arbitrary pure FCS. Finitely Correlated 
States are otherwise known as matrix product states or generalized valence-bond 
states. The bounds we derive become exact in the case where one considers the 
entanglement of a single spin with a half-infinite chain to the right (or the left) of 
it. Our bounds provide a proof of the recent conjecture by Benatti, Hiesmayr, and 
Narnhofer that their necessary condition for non-vanishing entanglement in terms of 
a single spin and the "memory" of the FCS, is also sufficient [6j. Our result also gen- 
eralizes the study of entanglement in the ground state of the AKLT model by Fan, 
Korepin, and Roychowdhury |13]. Furthermore, our result permits a more efficient 
calculation, numerically and in some cases even analytically, of the entanglement in 
arbitrary finitely correlated quantum spin chains. 

We continue the study of entanglement in the setting of ground states of Hamil- 
tonians with a spectral gap. In particular, for V a finite subset of Z'^, we let Hy 
denote a Hamiltonian on V with finite range, finite strength interactions and a unique 
ground state with a non-vanishing spectral gap. For a density matrix pA that de- 
scribes the finite-volume restriction to a region A of the unique ground state, we 



provide a detailed version of Hastings' proof in one dimension [23], that the entropy 
of Pa is bounded by a uniform constant C, where C depends only on the interac- 
tion strength, the spectral gap and the maximum among the dimensions of the state 
spaces associated with each site in V . Moreover, we provide a detailed generalization 
of the 1-dimensional construction of Hastings' approximation to the ground state in 
dimensions 2 and higher, which may prove useful for understanding the underlying 
structure of ground states of gapped Hamiltonians in higher dimensions. 

Finally, we turn our attention to the study of a conjecture central to Quantum 
Information Theory, the multiplicativity of the maximal output Schatten p-norm of 
quantum channels, for p > 1. In particular, we study the output 2- norm for a special 
class of quantum channels, the depolarized Werner-Holevo channels, and show that 
multiplicativity holds for a product of two identical channels in this class. Moreover, 
it is shown that the depolarized Werner-Holevo channels do not satisfy the entrywise 
positivity (EP) condition introduced by C. King and M.B. Ruskai, which suggests 
that the main result is non-trivial, since the EP condition has been shown to imply 
multiplicativity of the output 2-norm. 
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Chapter 1 



Introduction 



1.1 Informal Overview and Historical Motivation 

The newly developed fields of Quantum Information Theory (QIT) and Quantum 
Computation (QC) have cross-fertilized a variety of areas in long established fields of 
mathematics and physics, such as Operator Theory and Condensed Matter Physics. 
The mysterious resource of quantum entanglement, the stronger-than-classical non- 
local correlation between quantum particles, is at the center of many exciting new 
theories and applications, such as polynomial-time factoring and quantum teleporta- 
tion. The physical and mathematical framework behind many of these applications 
can be formulated in the language of one dimensional arrays of quantum particles 
(spins), known as quantum spin chains. These systems are completely described 
by their associated Hamiltonian, a mathematical object which encodes the interac- 
tions between the particles. More importantly, the Hamiltonian contains information 
about the possible states of the system it describes and the energy level associated 
with each such state. Of special importance is the state of the system with the lowest 
energy, often called the ground-state of the system. The minimum amount of energy 
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required to excite the system from its ground-state is called the spectral gap. 

Not long ago, quantum spin chains were primarily the object of study in Con- 
densed Matter Physics. Since the advent of QIT and QC, there has been a renewed 
interest in the properties of these seemingly simple quantum objects. The dynamics 
that govern spin chains can be applied to quantum states to transform them to de- 
sired target states, thus achieving quantum computation. The existence of a spectral 
gap in the Hamiltonian describing the dynamics of the spin chain is key to performing 
quantum computation reliably. The gap not only acts as a safeguard against exter- 
nal perturbations that might otherwise derail the computation, but it also determines 
the time it takes computational protocols such as adiabatic quantum evolution [2] 
to solve problems which are intractable for classical computers. Hence, the study 
of the spectral gap and its implications is central to efficient and reliable quantum 
computation. 

Recent progress in the study of the velocity with which interactions spread be- 
tween spins on a spin chain [IHl |21] has been combined with techniques used for split- 
ting Hamiltonians with a spectral gap into "frustration free" components |l3l [26] , to 
produce some spectacular results in Quantum Information Theory, with applications 
to Quantum Computation. A question central to the study of the complexity of sim- 
ulating quantum systems with classical resources [591 E] is the conjecture that the 
entropy of a bounded region of spins in the ground state of a gapped Hamiltonian 
grows proportionally to the surface area of the region and not the bulk (volume). 
This conjectured entropy scaling has attracted much attention lately as it may both 
help explain and optimize the running time of algorithms [SS], such as the Density 
Matrix Renormalization Group (DMRG) algorithm [63], that exploit the underlying 
structure of "frustration free" ground states, such as the one describing the ground 
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state of the well-known AKLT model |2T1|1], to calculate efficiently certain properties 
of a given Hamiltonian. The reason why such algorithms work so well at calculating 
properties of quantum systems is straightforward: The system, usually, has a lim- 
ited amount of entanglement between its components, hence it can be approximated 
efficiently by a series of ever-refined classical-like components, known as Finitely Cor- 
related States [ISllTS], or more recently, as Matrix Product States |60l[5]. This special 
class of states will continue to attract much attention, as we strive to understand bet- 
ter the conditions under which entanglement persists in the ground states of quantum 
systems we wish to use in practice. 

1.2 Summary of the Main Results 

In the chapters that follow, we present research that has either been published (Chap- 
ter [2] in IH] and Chapter |4] in [12]), or, is being prepared for publication (Chapter 

In particular, in Chapter |2] we present a proof of a conjecture posed by Benatti 
et al. in |6J. In their paper, using the structure of Finitely Correlated Pure States, 
introduced by Fannes et al. in [I5l [16], the authors construct a "dual" state F(p) 
that encodes a large part of the information found in the original state of the spin- 
chain. Using a brilliantly simple argument, they proceed to show that in order for 
entanglement to exist between a spin at position 1 of the spin-chain and a subset 
of spins in positions [p,n] {p > 1) in the original state, it is a necessary condition 
that the "dual" state F(p) has non-zero entanglement. Using numerical evidence, 
they conjecture then that this condition is also sufficient for entanglement between 
subsets of spins on the translation invariant state. In this chapter, we resolve this 
conjecture by proving the stronger claim that the entanglement of formation between 
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the first spin and the spins at sites [2, n] converges exponentially fast in n to the 
entanglement of formation of F(p). 

In the first section of Chapter |3| we turn our focus to the Area Law for 1- 
dimensional ground states of gapped Hamihonians and provide a detailed proof of 
a ground breaking result by Hastings [23], the proof of the conjectured area law for 



1-dimensional Hamiltonians with a spectral gap. In section 3^, we leave the realm of 
Quantum Information Theory and enter the field of Quantum Statistical Mechanics 
in order to construct a higher-dimensional generalization of Hastings' approximation 
for the ground state of a gapped Hamiltonian, which is central to the proof of the 
area law in one dimension. It is our hope that a clear presentation of the techniques 
and ideas involved in the construction of such an approximation, will shed more light 
on the underlying entanglement of the ground states of Hamiltonians that describe 
quantum systems with a spectral gap [131 EZl EH], leading to a better understanding 
of the scaling of entropy in 2 and higher dimensions. 

The results in Chapter |4] came about after a conversation with M. B. Ruskai at 
the 2006 International Congress of Mathematical Physics in Rio de Janeiro, where 
she posed a question related to the "multiplicativity conjecture" for a certain class 
of quantum channels. In mathematical terms, a quantum channel is a completely 
positive, trace-preserving map between the algebras of bounded operators on Hilbert 
spaces, and "multiplicativity" refers to the maximal output (Schatten) p-norm equal- 
ity 11$ ® ^E'llp = ||$||p ||\E'||p, where $, \E' denote quantum channels and ||fi||p = 
supp ||i7(p)||p, p being the input state. The study of this conjecture and its im- 
plications are at the center of Quantum Information Theory, as other important 
"additivity" conjectures involving the capacity of quantum channels and the role of 
entanglement in transmitting quantum information may be proven [561 119] by show- 
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ing the following sufficient condition: For each quantum channel $, there exists a 
p($) > 1 such that 11$ ® $||p„ = ||$||p„ for 1 < p„, < p($) and p„ -^ 1. Focusing on 
the special case of the maximal output 2-norm (a.k.a. Hilbert-Schmidt norm), I was 
able to prove multiplicativity for a class of channels given by $(p) = Ap+(1 — A) ^~^ , 



for all d > 2, A G [0, 1] [42j. The channels 1 studied are a generalization of the well- 
known Werner-Holevo channels W{p) = ^~^ , which gave the first counterexamples 
to multiplicativity (for p > 4.79). It was recently shown that for p y^ 2 there exist 
quantum channels that yield counterexamples to the conjecture of multiplicativity 
[6^ [27] . To the best of my knowledge, there is no published counterexample for 
p = 2, but the focus has shifted to studying multiplicativity in the region of p < 1. 
In any case, like with so many other questions in Quantum Information Theory, the 
study of multiplicativity and entanglement in general, promises to be both fruitful 
for various areas of mathematics, physics and computer science and truly exciting in 
the development of both the theory and future applications! 



Chapter 2 



Entanglement in Finitely 
Correlated Spin States 



2.1 Introduction 

Entanglement properties of quantum spin-chains have recently attracted attention 
from researchers in quantum information theory and condensed matter physics. From 
the perspective of quantum information theory, the distribution of entanglement over 
long ranges via local operations on a spin-chain [HI [121 IS] has obvious applications 
to teleportation-based models of quantum computation [S2l CHI EO]- Moreover, it 
has recently been shown that entanglement in finitely correlated chains [15j can be 
used to achieve universal quantum computation ^60j and provide a computational 
tool for adiabatic quantum computation [5]. On the other hand, the scaling behavior 
of quantum correlations in infinite spin-chains is intimately related to their critical 
behavior (recent work has established a general mathematical framework for studying 
entanglement in infinite quantum spin-chains [33].) 
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Finitely Correlated States (FCS) are a generalization of the so-called Valence Bond 
Solid (VBS) states, which arise as the exact ground states of a considerable variety 
of quantum spin Hamiltonians |^ [T]. Interestingly, FCS also provide approximate 
ground states of any quasi-onedimensional spin system with finite-range interactions 
[17] . In particular. Density Matrix Renormalization Group calculations produce nu- 
merical approximations of the ground state of spin chains that can be regarded as 

FCS isn]. 

Motivated by the potential applications of distributed entanglement in finitely 
correlated chains, Benatti, et al. in [6], give a necessary condition for entanglement 
between a spin and a subset of other spins; namely, that the entanglement between 
a spin and, what one may think of as [37j, the "memory" of the finitely correlated 
state must be non-zero. They, furthermore, conjecture that the same condition is 
sufRcient, in the sense that it implies entanglement between a spin and a subset 
of other spins. We present here a proof of that conjecture by showing that the 
entanglement between a spin and its neighbors converges exponentially fast (in the 
number of neighboring spins) to the entanglement between a spin and the "memory" 
of the finitely correlated state. Moreover, we show that entanglement between distant 
spins vanishes exponentially fast in the length of their separation. 

Since finitely correlated states provide the exact ground states for generalized 
valence-bond solid models [33], our result generalizes the calculation of entangle- 
ment [T3] for the AKLT model [I]. 

More importantly, our result implies a simple and computationally efficient way for 
detecting distributed entanglement in finitely correlated states. Namely, the Positive 
Partial Transpose (PPT or Peres-Horodecki) criterion [CT ISO] can be applied to the 
state describing the interactions of a spin with the "memory" of the finitely correlated 
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state, to detect entanglement between a spin and a subset of other spins. 

2.2 The setup and main result 

We will work with translation invariant pure FCS [15] on the infinite one-dimensional 
lattice. For each i G Z, the spin at site i of the chain will be described by the algebra 
A of d X d complex matrices. The observables of the spins in an interval, [m, n], 
are given by the tensor product A[m,n] = ®'j=mi^)j- The algebra Az describing the 
infinite chain arises as a suitable limit of the local tensor-product algebras A[-n,n] '■= 
®"=_„(^)j. Any translation invariant state u over Az is completely determined by 
a set of density matrices P[i,„], n > 1, which describe the state of n consecutive 
spins. In the case of a pure FCS, as was shown in \16l, these density matrices can be 
constructed as follows: 

The memory, B, of a FCS is represented by the algebra oi b xb complex matrices. 
Let K : A ® B ^-^ B he a completely positive unital map of the form K{A ^ B) = 
V{A ® B)V^, where V" : C^ ® C^ ^ C^ is a linear map such that VV^ = 1b- We 
define the completely positive map K : B ^-^ B,hj E,{B) = E(lL_4(8)i?). The condition 
on V implies that E is unital: E(lLg) = Ug. In y^ it is proved that for any pure 
translation invariant FCS, it is always possible to choose B and V such that there is 
a unique, non-singular bxb density matrix, p, with the property Trp]E(i?) =Tt pB, 
for all B eB. 

We introduce the density matrix pa(^b associated with the state encoding the 
interaction between the spin at site 1 and the "memory" of the FCS: 

Tr^0H {pA<^B A0B^= Tre (p E{A ® 5)) . 
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Using the cychcity of trace we also have pA(g)B = V^pV. 

We are now ready to define the density matrices P[i,n,] recursively, by the following 
identity: 

Tr(p[i,„]Ai® ^2®- ■ ■» A„) = Ttb{pE{Ai0E{A20-- ■®E(A„_i®E(A„® 1^)) ■ • ■ ))). 

(2.2.1) 
From the above definition and the cyclicity of the trace we get the equivalent definition 

P[i,n] = TrB(K!p^®BK), 

where K = (1^ ® V^)(l^®2 ®V)--- (U^^n-i V). 

An important property, intimately related to the exponential decay of correlations 
in a pure FCS is that the peripheral spectrum of E is trivial; that is, Ug is the only 
eigenvector of E with eigenvalue of modulus 1 [16j. This implies that the iterates of 
E converge exponentially fast to E°° given by E°°(i?) = lim„_,ooE"(i?) = Ti(pB)l.s. 
More precisely, for any A such that |Aj| < A < 1, for all eigenvalues Aj of E different 
from 1, there exists a constant c such that for all n > 1: 

||E"-E°°|| <cA'", (2.2.2) 

where the norm is the oo-norm on B considered as a Banach space with the 1-norm. 

Our object of study is the entanglement of formation, EoF [7j. The EoF is defined 

for states of composite systems with a tensor product algebra of observables ^"1(8)^2. 

Definition 2.2.1 (Entanglement of Formation). The entanglement of formation of 
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a bipartite state over Xi ® X2 with associated density matrix au is given by: 



E[x,,X2]i(^i2) = mi^piSlTixMl 



12JJ^ 

where S{p) = —Tip log p is the von Neumann entropy and the infimum of the average 
entropy is taken over all convex decompositions ai2 = J2i Pi ^12 i^^o pure states. 

Whenever Xi is finite dimensional, as will be the case for us, the infimum can be 
replaced by a minimum in the above definition, i.e., there is an optimal decomposition, 
{pi, crl2}, where the infimum is attained (see [10] for details). We call {|0j)} an 
ensemble for the density matrix a whenever the latter can be decomposed as a = 
Yli\4'i){4'i\- There are an infinite number of ensembles corresponding to a given 
density matrix. The following lemma provides us with a complete classification: 

Lemma 2.2.2 (Isometric Freedom in Ensembles, [MIES])- Let {|ej)}f^]^ be the en- 
semble corresponding to the eigen- decomposition of the density matrix a, where d = 
rank{a). Then, {|^i)}^i is an ensemble for a if and only if there exists an isometry 
f/ : C^ ^ C™ such that 

d 

\^i) = X^ ^'.i l^i)' 1 < ^ < "^• 

The above lemma implies that any two ensembles for the same density matrix, 
{\il)j)}j^i., and {|0j)}j=\, are similarly related via a partial isometry W : C*^^ 1-^ C^^, 
Our main result is the following theorem: 

Theorem 2.2.3. For any pure translation invariant FCS we have 

< E^j^^t3]{.pA®B) - ^[A^»"-i](P[i,n]) < e(?^), (2.2.3) 
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where e{n) decays exponentially fast in n. 

2.3 Proof of the Theorem 

The lower bound is proven in [6]. For the sake of completeness, we include here the 
following proof. 

The definition of p[i^n] implies that every decomposition of p_A(S)B into pure states 
induces a decomposition of P[i,n]- Moreover, the restrictions to the spin at site 1 of 
the i-th state in the corresponding decompositions of pa(^b and p[i_„] are equal. To 
see this, note that since the operators Vn leave the first spin invariant, the cyclicity 
of the trace implies 



Tr^®n-i(p[i_„]) — Tr^®n-icge(V„p^jggK) — Tri3(p; 



A(^B)y 



where we have used VnV^ = 1^ ® Ig- It follows that for each decomposition of pa(^b 
there is a corresponding decomposition of p[i,„] with equal average entropy. Since the 
average entropy of P[i,n] is minimized over a (possibly) larger set of decompositions, 
the lower bound follows. 

We now focus on the upper bound. We start with the following decompositions 
of Pa«iB and p[i^n] into (unnormalized) pure states: 

6 

PA^B = 1]V^^|X.)(X«I^ (2.3.1) 

b 

P[i,.] = J2^n,^^\^^)^^^\^^-d, (2-3-2) 

where {|Xj)}i=i is the eigen-ensemble of p and Gnj = Ki(lL_4«n (g) |Xj)/IIXill)- The 
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term in parenthesis in the expression for Gn.j comes from the Kraus operators in the 
decomposition of the completely positive map Tr^. 



By the observation following Lemma 2.2.2 we have that the (unnormalized) states 



|$") in the optimal decomposition of P[i,n] are given by: 



l'^r) = Ef^'.foA.^^lx.)' ^<1<L, (2.3.3) 

for some partial isometry f/ : C^ ^-^ C^, whose dependence on n we suppress. 
Moreover, it is easy to check that pa(^b has a decomposition into J2i'^i^ with ai = 
E,I^Kj))(^/(j)|and 

b 

I^Kj)) = 5Z«a..) VPiV^X^), 1<1<L (2.3.4) 

To calculate the EoF we need the restrictions of {|$")($"|} and {a/} to A: 

0r = Tr^«"-i(|$r)('^ri), ai = Tre(^0- (2.3.5) 



Define the density matrices 0" = ^"/a" and ai = <Ji/Pi, where a" ; 

Tr(0r), A^ll^/lli =Tr(aO. 

From the definition of the EoF and the optimality of {0"}^]^ we get: 



^nii 



L 



E[A,t3]{PA(s>B) - E[A,A(s^-^{pii,n]) < ^e^^)' (2.3.6) 



1=1 



where ei{n) = piS^ai) - ^'^(^r)- 

It remains to show that Yli=i ^i{^) is exponentially small. We estimate each term 
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in the sum as: 



Hn)\ < (3i\S{ai) - 5(0r)| + lA - an ^ogd, 



(2.3.7) 



since rank(</)") < d. 

To bound \S{ai) — S{(f)]^)\ we use Fannes' inequahty for the continuity of the von 
Neumann entropy [18]: 



\Siai) - S(0r)| < (logrf + 2)\\ai - 0r||i + r/(||a, - 0rili), (2.3.8) 



where ri{x) 
we have: 



— xlogx and log is the natural logarithm. By the triangle inequality 



iA-ari = iii^dii-ii0riiii<ii^^-0riii- 



(2.3.9) 



Another application of the triangle inequality gives: 



\(^i-0i\\i 



< 



A 



which simplifies, with the use of (2.3.9), to the following inequahty 



Wi-<P?\\i<2 






(2.3.10) 



Combining equations (2.3.7)-(2.3.10) and setting 



7^"^iia,-0riii/A, 



(2.3.11) 



we get the following bound for e;(n): 



|Q(n)|<A[(logrf3 + 4)r," + r/(2rr)]. 



(2.3.12) 
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where we have assumed that 2r" < 1/e, to assure ri{x) is increasing. 

To complete the proof, we show that r" is exponentially small for large n. Since 
each Gnj leaves the spin at site 1 invariant, the cyclicity of the trace yields: 

b 
'^r= E Uli^'f)U^^,)^^B{V^X^){X^'\VG^,fGl,), 

But Gn,fGl^ = l^®E"-i(|x,')(x,|)/(||x,'|lllx,ll)- Substituting E- for E"-i we get: 

b 

where X,,, = V^x^){x^'\V and F,,,. = U^ ® [E- - E"-i](|x/)(x,|)/(||x/|| ||x,||). 

Like all trace preserving quantum operations, the partial trace is contractive with 
respect to the 1-norm. Hence, an application of the triangle inequality for the 1-norm 
gives: 

b 

||5-/-0r||i< 2^ |f^it(*'i')ll^''fo')lll^«.^'llill^i'i'lli 
It is not hard to see that 

II \^ II lU, IIIU, II ||"V^ II ^ Il11?("-— 1) wooll 



and hence 

b b 



riii<(EEif^'-fe-)iiix.ii) iiE("-^)-E° 



1^/ 
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Since ^- . ^ |t^i,(ij)| llXill — A? two apphcations of Cauchy-Schwarz give: 



m-<Pi 



nil <6^AiiE^""^^-E°°i 



(2.3.13) 



Finally, combining (2.2.2) with (2.3.13), equation (2.3.11) becomes 



rr < CiA", ci = cbyx. 



(2.3.14) 



To conclude the proof, we note that since the bound for r" is independent of I, 



summing over / in equation (2.3.12) yields: 



J2 K\ < (log^' + 4)ciA" + r/(2ciA"). 



1=1 



It is clear that for A' > A there exists a constant C2 such that 



r/(2ciA") < C2{X'y 



The only condition on n was imposed in equation (2.3.12) were we assumed that 



2r" < -. Using equation (2.3.14) we see that there is an no such that the above 



condition is satisfied for all n > uq. The previous observations imply that for all A' 
with A < A' < 1, there is a constant C3 such that: 



e[n) = C3 



(A')"> J]|eri, for all 



n. 



1=1 



Finally, equation (2.3.6) implies that: 



E[A,B(g,B]iPA(E>B) - E[A,A'^"-^iP[l,n]) - ^l^)' 
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and this completes the proof of the theorem. 

A natural question to ask at this point is the following: How does the entanglement 
between the spin at site 1 and spins at sites [p,n], p > 2 behave as p becomes large? 
Since the state Pi,[p,n] factorizes into pi ® P[p^n] as p — > oo [28], we expect that the 
bulk of the entanglement is concentrated near site 1. The following theorem confirms 
this: 

Theorem 2.3.1. For any pure translation invariant FCS andn > p >2, the following 
bound holds: 

E[A,A'^"-p+^]ipi,[p,n]) ^^ lnd{n-p + l)e{p), (2.3.15) 

where C is of order unity, d is the dimension of each spin and e{p) decays exponen- 
tially fast in p. 

Note that the above theorem implies exponential decay of the entanglement be- 
tween spin 1 and spins [p, n] as long as n — p does not grow exponentially in p. 

Proof. The main observation is that the trace distance between the states Pi,[p,n] 
and pi (g) pip^n] vanishes exponentially fast with p. This is a consequence of the 



exponential rate of convergence described in equation (2.2.2). To see this, note that 
IIPi, [p,n] - Pi ®P[p,n] 111 = Tr[(pi^[p^„] -pi (g)p[p^„]) P], where P is the projection onto the 
positive eigenvalues of pi,[p,n] — Pi ® P[p,n]- But, 



Tr(pi,[p,„] P) = J2 Tr(p[i,„]P,^®U[2,p_i]®i^,^®---®P,J 
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and applying (2.2.1) we get, 



TT{pi,lp,n] P)= Yl '^^^ {pHPh ® EP-2 ( E(P,^ ® ■ ■ • ® E(Pi„ 1b)) ■ ■ 

(2.3.16) 
where P = V ■ • ■ ■ _, Pi, ® -Pj„ ® ■ ■ ■ ® Pi,, is a decomposition of P into sim- 

.6 ' '1 (t-p )''p+l V5*T1 -L-L P "■ ■*■ 

pie tensor products over some fixed basis for each tensor. Furthermore, from the 



definition of E°^ and (2.2.1) we have 



Trpi®p[p,„,]P= J2 TrB(pE(Pi, ®E°°(E(P',^®---®E(P,„®Ub))-- 

^1 1^p?'^p+l V!^^^-'- 

(2.3.17) 



Combining (2.3.16) and (2.3.17) we get that ||pi,[p,n] — Pi ® P[p 



i,n]\\l 



^1 )'^P)^p4-l !•■ -i^n — J- 



J2 Tr^ (pE(P,, ® [W-^ - E°°) (E(P,^ ® 

< ||U^®(EP-2-E-)||^ |^TrB(pE(P,, ®(E(P,^« 

< cAP-2Tr(p[i,„_p+2]P) <cAP-' 



E(P,„®Ub))-- 

>§E(Pi„®lB))---) 

(2.3.1 



Now, note that for two density matrices p and a the (normalized) trace distance 
T{p,a) = ''^~°"''^ is an upper bound on the Bures distance D{p,a) = 2a/1 — P(p, cr) 



[9], where P(p, cr) = Tr y p^/^crp^/^ is the fidelity measure. In particular, the following 
bound holds |1H1 Ch. 9]: 

P>(p,(t)<2v/T(p,(t). 

Since -E[yi,y!i®n-p+i](pi®P[p,„]) = 0, a straightforward application of Nielsen's inequality 
for the continuity of the EoF |19] yields the desired result. D 
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2.4 Discussion 

Having established such a strong connection between the states p[i^n] and pa^Bi one 
can apply various entanglement criteria on p^^g to deduce entanglement properties 
of the spin chain. To start with, we note that for qubit chains with 2-dimensional 
memory algebra B, the entanglement of P[i,n] can be computed analytically (in the 
limit) by evaluating the concurrence [22] of Pa®b- For higher dimensions one can 
apply the PPT criterion to Pa^b to detect distributed entanglement in the finitely 
correlated state. Specifically, the main theorem in [31] implies that there can be 
no PPT bound entanglement in p^^g since rank(p^0g) = b < inax{d,b}. Hence, 
if the partial transpose of pa(^b is positive, then p^^g is separable. On the other 
hand, if the partial transpose of pa(^b is negative, then for n large enough p[i_„] 
becomes entangled. The amount of maximum entanglement in P[i,n] depends on 
the amount of entanglement found in pa®b- From this point of view, it would be 
very interesting to look at FCS that maximize entanglement of Pa®b- Moreover, 
understanding how entanglement of PA(^b varies with different CP maps E could lead 
to a better understanding of how phase transitions occur when we vary the parameters 
in the underlying Hamiltonian of the system. 

To conclude, we note that the conjecture of Benatti, et al. [6], follows as a corollary 



of Theorem 2.2.3 In particular, our result implies that a spin at site 1 of the chain 
is entangled with spins at sites [2,n] (for n large enough) if and only if p_a®b is 
entangled. Moreover, the entanglement of P[i,n] approaches the entanglement of pa®b 
exponentially fast. 
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Chapter 3 



Entanglement in the Ground State 
of Gapped Hamiltonians 



3.1 Introduction 

It is widely believed that the entanglement entropy of a region A in the ground state 
of a gapped quantum spin Hamiltonian with finite range interactions does not grow 
faster than the area of the boundary of A. Although this property in its general 
formulation is still a conjecture, this is called the Area Law for the entropy. The Area 
Law is of interest not only for theoretical reasons, but also because it has practical 
implications for the computational complexity of calculating the ground state with 
a desired level of accuracy [61j. The intuition behind the Area Law is simple. By 
the Exponential Clustering Theorem, a non-vanishing spectral gap implies a finite 
correlation length |l5l |22] and this puts an exponentially decaying bound on the 
entanglement of two spins as a function of the distance. From this, one may guess 
that only spins near the boundary contribute significantly to the total entanglement 
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with the exterior of the region and hence, an area law should hold. The relationship 
between the correlations and entanglement of a region with its exterior, however, is 
not sufficiently well understood to lead to a proof based directly on this intuition. 
If one assumes a decay property of the mutual information instead of correlations, 
one can indeed prove an area law [B5j. The Area Law itself has only been proven 
in one dimension by Hastings [23]. In more than one dimension, for special systems 
with valence bond (i.e. matrix product) ground states, the area law is easy to derive, 
but there is no general result. Here, we review in detail the one dimensional result 
by Hastings and provide a higher dimensional generalization of his ground state 
approximation. Although the Area Law bound we discuss below applies only to 
1— dimensional systems, we describe the result in a more general setup in case this 
leads to future insight on how to treat the Area Law question in higher dimensions. 

3.2 Setup and Main Results 

We will consider a finite system of spins located at sites in V, a finite subset of Z*^. At 
each X G V, we have a finite-dimensional Hilbert space of dimension Ux and assume 
that rix < A^, Vx G V, for some constant A^. Let Ax be the algebra of observables 
associated with X G V. We consider the dynamics tj generated by the Hamiltonian 



Hy=Y, HX). 



xcv 



The main assumptions are as follows: The interactions $(X) are uniformly bounded 
and finite range, and, for convenience, we will take pair interactions with range 1. 
So, for all X gV, $(X) = $(X)* G Ax, \\^x\\ < J, for some constant J > 0, and 
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$(X) = if diam(X) > 1. Here, for X ^V, 

diam(X) = max{(i(a;, y) | x, y G X}. 

We will also assume that the Hamiltonian Hy has a unique, normalized ground 
state, which we will denote by l^&o), and a spectral gap 7 > to the first excited 
state. 

For a set F C V , the boundary of F, denoted by BY , is 

dY = {xeY \^y (^V\Y, with d{x, y) < 1}. 

Moreover, for ^ > 1, we define the following sets: 

ly = iy{C) = {x G r I Vy G (9F, d{x, y)>l}. (3.2.1) 

The set ly corresponds to the ^-interior of Y and it will be empty if £ > diam(y). 

By = Brie) = {xeV \3yedY, d{x, y) < i} . (3.2.2) 

The set By corresponds to the ^-boundary of Y. 

Ey = Ey{i) = {xeV\Y\yyedX, d{x, y) > i} . (3.2.3) 

The set Ey corresponds to the ^-exterior of Y. Note that ly, By and Ey are disjoint 
and moreover, V = ly U By U Ey, as can be seen from Figure |3.2.1[ 

Finally, define Bi = Bj{i) = Y n By{i) and Be = Be{(.) = {V \Y) f] By{i) 
and set Dj = HxeBj ^x and De = ^x^Be ^x- We will assume from now on that the 
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Figure 3.2.1: The different regions around tlie boundary of A. 



volume of tlie interior of tlie ^-boundary satisfies 



\Bi\ <re\dY\, 



for some constant r > 1. For convex Y, tliis inequality liolds witli r = 1. Tlie regions 
we will be working with will be convex. 

We denote by Pq the orthogonal projection onto |\E'o) s-nd by p^ the density matrix 
describing the ground state restricted to the region A: pA = Try\^Po- 

We are now ready to state the main result of 



Theorem 3.2.1 (Area Law). Let Hy denote a 1- dimensional Hamiltonian on an 
interval [1, \V\\, with the properties described above. For any 1 < M < \V\, let 
A = [1,M]. Then, the following bound, independent of M applies to the entropy of 
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the region A: 

S{pa) <C^lnC{-f,d,J)lnN2^^^''^. (3.2.4) 

where C is of order unity and the constants ^ and C(7, d, J) are defined in Theorem 



We will provide the complete proof of Theorem 3.2.1 assuming Theorem 3.2.2 



below, which is a generalization to higher dimensions of a result of Hastings 



The proof of Theorem |3.2.2| proceeds along the same lines as in the one-dimensional 
case [221 1251 [21] , but it is a bit long and technical and therefore we present its proof 
in section [3^ 

Theorem 3.2.2 (Ground state approximation). There exists .^ > 0, such that for 
any sufficiently large i > cq d^^'^ and Y G V , there exist two orthogonal projections 
Py G Ay , Pv\Y ^ -A-vXY o^nd an operator Pb G Asyii) ^^^h \\Pb\\ < 1; such that 

\\PbPyPv\y-P4 <C(7,rf,J)|ar|V^/« (3.2.5) 

where Pq is the projection onto the (unique) ground state, cq is of order unity and 
C{''y,d,J), C, are explicit in terms of J, 7 and d, the dimensionality ofV. 



The operator Pb in (3.2.5) is responsible for all correlations in the ground state 
approximation between the region Y and its exterior. Its support is concentrated 
along the boundary of Y . This is reminiscent of the structure of matrix product states 
[H Uni |60] . The finite extent of correlations and entanglement across any boundary 
is essentially a consequence of the non-vanishing spectral gap and the existence of 
a finite Lieb- Robinson velocity [5S1 SSI [221 US]- The problem of calculating such an 
approximation of the ground state is a related but separate question we intend to 
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turn to at a later occasion. A numerical algorithm for this problem is discussed in 
[26] . It is known that in general this is an NP-hard problem |13j . 

3.3 Proof of Main Theorem 



We present now the proof of Theorem 3.2.1 



Proof. First, we introduce the length m* = ^i^^fg^ and then set Y = Y(m) = 
lA{m), m < m*. The first step is to use the well known subadditivity of the Von 
Neumann entropy, which follows from the non-negativity of quantum relative entropy 
S{p\\o-) = Tr(plnp — plncr), to get the following bound for S{pa)'- 

S{pa\\py ® pa\y) = S{py) + S{pa\y) - S{pa) > (3.3.1) 

The above bound, combined with the fact that 

S{pa\y) < |A \ r| IniV < m*\dA\ \nN = S{pa)/^, 

implies that 

4 
S{pa) < - S{pY{m)), m<m* (3.3.2) 

We can now focus on bounding the entropy of the region Y{ra). We will do this 
by looking at the rate of decay of the Schmidt coefficients in the Schmidt decom- 
position of |\E'o) along the boundary of Y{m). Since we will be using the Schmidt 
decomposition of the ground state \^o)^ let us introduce it here as 



l^o) = Y^ Vao(a)|^y,o(a)) ® |^y\y,o(a)), 
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where J2a'^o(^) ~ ^ ^^^ {I'^Yfiict))}, {|^v\y,o(Q^))} are orthonormal sets supported 
on Y and V\Y, respectively. We order the Schmidt coefficients of |\E'o) in decreasing 
order such that if a < /5 then o-Q^a) > o"o(/3). Moreover, 

a 

and 

PV\Y = 5^ao(/5)|^y\y,o(/5))(^y\y,o(/3)|. 

Having fixed m* and y(m), we are interested in the overlap of the density matrix 
Py ® Pv\Y with Pq, since we will be treating the cases when the overlap is small and 
when it is large, separately. We define the overlap to be 

P = P{m)=TT{PopY(»Pv\Y)- (3.3.3) 



We use the approximation operator PbPyPv\y of Theorem 3.2.2 in order to relate 
the overlap P with the approximation error e(£), over which we have some control. 
Thus, we choose the following approximation to the ground state: 

p{i) = PbPyPv\yPy ® Pv\yPyPv\yPI (3.3.4) 

which, as we will see shortly, has Schmidt rank bounded by an exponential in m and 
\dY\, but independent of the volume of Y . More concretely, we will show that p(£) 
has a decomposition into pure states, each with Schmidt decomposition along the 
boundary dY with Schmidt rank at most A^^^l^^l. To see this, first note that pure 
states in the decomposition of PyPv\yPy ® Pv\yPyPv\y are product states along 
dY. Using the spectral decompositions of py and pv\Y introduced earlier, we may 
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focus our attention to product states of the form 



|^y(7, 6)) = Py|^y,o(7)) ® Py\y |^V\y,0 (5)) • 



(3.3.5) 



We study now how the action of Pb on |\l/y(7, S)) affects its Schmidt rank. 

Since Pb is an operator acting non-trivially only on sites in a subset of Byii), we 
have the following general decomposition: 



Dr 



Pb= Y^ 1/y ® E{a, P) ® G{a, /?) ® l^^ 

a,(3=l 



(3.3.6) 



where the De x De matrices G{a, j3) act on sites in Be and the matrix units E{a, (3), 
which act non-trivially on Bj, form an orthonormal basis for Dj x Dj matrices. 
Moreover, Dj < iVl^^l < A^^l^^l. 

To bound the Schmidt rank of PB|\E'y(7, 5)) we trace over sites in Y and study 
the rank of the operator 



Try (^PB\^Yil,5)){^Yil,5)\Pl) =Y, 



0=1 



Dr 



J2 c{p,p')\F^m{UP') 



./3,/3'=l 



(3.3.7) 



with \F^{P)) = G{a,P)0lEyPv\Y\'^v\Y,oi5)) and c(/?,/?') = {^yMIPy^Iy 



i?(/?', /?)Py I \E'y,o (7)), coming from (3.3.6) and the definition of |\E'y(7, 5)). 



Clearly, as a sum of Dj matrices each with rank at most Dj, the above operator 

has rank bounded by Dj and, hence, the Schmidt rank of PB|^y(7,'^)) is bounded 
aboveby Ar2£|9y|_ 

Now that we have a good grasp on the entanglement overhead produced by the 
approximation operator PbPyPv\y, we return to the question of how the overlap 



P relates to the approximation error e(£). Remembering from (3.3.3) that P 
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Tr(-PoPy ® Pv\y), we have the following overlap estimate between p(£) and Pq: 



v/Tr(PoPW) = \\PoPBPYPv\Y^/p^®,/p^h 

> \\PoVp^ ^ ^/pV^h - \\Po{Po - PbPyPv\y)\\\\Vp^ ® ^/pr^h 

> VP-e{i), 

where ||X||2 = vTrXXt and ||-|| is the sup-norm. Hence, for all £ such that e^(£) < P, 
we have 

Tr(PopW)> [^/P-e{i)y. (3.3.8) 

Moreover, we have that 



Tr[(l - Po)p(£)] = \\{1 - Po)iPBPYPv\Y - Po)Vp^ ^ y/fH^M 

\\Y 



< \\1 - Pof\\PBPYPv\Y - Pof < e\i) (3.3.9) 



Finally, upon normalization the overlap becomes: 

Tr(Pop(f)) _ 1 1 e\i) 

Tr(p(£)) " ^ " 1+ Tr(Popffl) ^ 1 - - u^i;^ ^ ^ - 2^^ (3.3.10) 



We will now use (3.3.10) for the case when P > 2e^(m*). The case P < 2e^(?7i*) will 
be treated later on. 

Case I: P > 2e'^{m*). 

In this case, the next step is to relate the above overlap to the Schmidt coefficients 
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of |\I/o)- More specifically, we will now show that for i > 0: 

To prove this, first note we already showed that p{i) is a convex combination of pure 
states PB|\t'Y(7, 5)), each with Schmidt rank bounded above by A^^^l'^^l. Let |^o(^)) 
be the (not necessarily unique) pure state in the aforementioned decomposition of 
p{i) satisfying 

Tr(Po|^oW)(^oWI) > ^^j^ (3-3.12) 

and introduce its Schmidt decomposition as 

s s 

J^V^l'^vm ® \<^v\Y{f3)), with J]r^ = 1 and s < iV^^I^^I, (3.3.13) 

/3=1 13=1 

as we have already demonstrated. For notational convenience, let My(Q;,/5) = 
|($y(/3)|^y,o(a))| and Mv\Y{a,P) = | ($y\y(/9)|^v'\y,o(«))| • Note that since each 
of {|$y(/?))}, {|$y\y(/?))}, {|^y,o(a))} and {|^y\y,o(a))} is an orthonormal set, 
Bessel's inequality implies that 



^My(a,/5)2 < 1 and ^Mv\Y{a,pf < 1, 

/3=1 a 

as well as 

^My(a,/3)=^ <1 ^ ^My(a,/3)2<s. 

Then, an application of the triangle inequality followed by Cauchy-Schwarz gives 
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the following upper bound for Tr(Po|*oW)(^oWI) = |(^oWl^o)P 



K^oWl^o)!' < iY,VMa)VW)MY{a,(3)Mv\Y{a,(3) 

\a,f3 

< (5^ao(«)My(a,/3)M (^r(/3)Mv\y(a,/5)^ 

V a,/3 / \ a,f3 



The last inequality follows from Schur convexity of f{[p{a)]) = Xla^o(tt)p(tt) and 
the observation that the vector [1, 1, . . . , 1, 0, . . . , 0], with at most s ones, majorizes 
ZpMY{l,(3r,EpMy{2,(3f,..: . 

To see that f{[p{a)]) is Schur convex, note that if we set Sp{a) = Ylk=iPi^) 
and A{a,P) = ao{a) — <Jo{P) then the condition that [p{a)] majorizes [q{a)] {p ^ q) 
becomes p ^ q 4^ Sp{a) > Sq{a), Va. Moreover, 



/( w«)}) - /({?(«)}) = E ^("' " + 1) (^p(") - ^^(")) 



which is non- negative since we have arranged the cro(a) in decreasing order so that 
A{a,a + 1) > 0, Va. 



Now that we have demonstrated (3.3.11), we may use it in combination with 



(3.3.10) to show that S{py) satisfies an area law in Case I. We begin by setting i' 



to be the smallest integer i such that 2e^(£)/P < 1. Using (3.3.10) and (3.3.11), we 
have for i > i' 

J2 ^o(«) < exp[-2(£ - i')/^]. (3.3.14) 
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We now maximize the entropy of py = Try\yPo5 given by 



Sipy) = -^ o-o(tt) ln(o-o(a)) 

a=l 



subject to the constraint (3.3.14). Following the notation of Lemma 3.3.2, set s„ 



j^2(i'+n)\dY\^ n > 1 and c = exp[-2/e]. Then, R = AT^I^^I and si = N^i^'+^Y] ^nd 



Lemma 3.3.2 implies that S{py) is bounded above by: 



{^ + 2 + 2i') In N\dY\ + [ ^ + I)ln2, 



(3.3.15) 



where we used the inequalities 2c/(l — c) < ^, 1/(1 — c) < ,^/2 + l and H2{l — c) < ln2. 
Using the definition of i' we have that 



2,- < ^ ,„ I ?5!(^4^)W V 2 ^ ^ ,„ ^2^ I + 2,„. + 2. 



P 



P 



Since we are considering Case I, we have furthermore that 2£' < 2m* + 2. Hence, 



using (3.3.15) we get 



SipY)<Ci\nN\dY\ +SipA)/2 + C2 



(3.3.16) 



with Ci = ^ + 4 and C2=(| + l)ln2. Using (3.3.2) and rearranging terms gives 



SipA)<^Ci\nN\dY\ +4:C2<4:{Ci\nN + C2)\dA\, (3.3.17) 



which is an area law bound. 
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It remains to treat the following case: 



2/™*^ 



Case II: P < 2e\m _ 

The main idea here is to show that in this case, there is a non-trivial lower bound on 
the relative entropy 

Sr{£) = S {pB(e)\\pBiie) ® pBE(e)) = S{pBi{i)) + '5(pbe(£)) - S{pb{i)), (3.3.18) 

remembering that B{1) = By{€) is the £— boundary of Y{m) = I Aim), for m < m*, 
and that Bj{i), BE{i) represent the interior and exterior of B{i) along dY, respec- 
tively. In particular, we will show that for any fixed m and for i sufficiently large, 

SR{i) > -(1 - e{i)f ln(2P + 4e(£)) + (1 - e{i)f ln(l - e{i)y - ln2, (3.3.19) 

which, in the case that e{m*) < 1/2 and P < 2e^(m*), implies 

Snii) > -{l-e{£))Hn{Ae\m*)+4e{e)) + {l-e{£)f\n{l-e{£)f-\n2 

> -(1 - e{£)y\n{2e{m*)+4e{£)) + (1 - e(£))Mn(l - e{£)f - In 2 

> -(1 - e{£)y \n{Qe{£)) + (1 - e{£)f ln(l - e{£)f - ln2 

> -lne(£)-lnl2, £<m*, (3.3.20) 

where we used (l-e(£))2ln(l-e(£))2 + (l-(l-e(f))2) In 2e(£) > -H2{{1 - e{£))^) > 
— In 2 in the final inequality. Note that if e{m*) > 1/2, then we have the following 
upper bound on m*: 

m*<ein(2C(7,rf,J)|9rn, 
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which imphes an area law bound on S{pa) with a logarithmic correction: 

SipA) < (4ein(2C(7,rf,J)) + 8^1n|M|)lnAr|5v4|, 

where we have used the definition of m* and the simple inequality | dY \ < \ dA \ . 



We now go back to proving (3.3.19), which, combined with (3.3.18) and (3.3.20), 



ultimately implies for m < m* and ii < i < m*: 

S{Pbm) + 5(pB,(,)) - S{pBi^)) > ^ - ln(12C(7,rf, J) \dY{m)\') = ^, (3.3.21) 

where ii is the smallest £ for which e{£) < 1 and £2 = ^ln(12C(7, (i, J) |9F(m)p) < 
^lnl2 + £i. 

We begin by restating the following important result about the relative entropy 



Lemma 3.3.1 (Monotonicity of Relative Entropy). Let Ti and /C he finite di- 
mensional Hilbert spaces, and let p, a he two density matrices on Ti. Then, for any 
trace preserving, 2-positive map $ : BiH) -^ B{1C), the following inequality holds for 
the relative entropy: 

S{p\\a)>S{<^{p)ma)) (3.3.22) 



In particular, for any operator M on Ti, with ||M|| < 1, setting 

Tr((l-|M|)p) J' 

we get a trace preserving, 2-positive map $ : BiTi) -^ Biil?). To verify that $ is 



^, , . Tr(|M|p) 

$(p) 



(3.3.23) 
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2-positive, we note that for any positive 2x2 block matrix A 



the determinant of the 4x4 matrix U2 ® $ (A) 



^ $(a) $(6) ^ 



a 



> 0, 



V 



$(6t) <l>(c) 



is given by the 



product of the determinants of the following two matrices: 



/ 



Uo^Tr 



U2®Tr 



/ 



VW\ 



V 

// 



\ ,/\M\ J 



a 



/ 



M 



\ 



v 



v 



'1 - \M\ 

Vl-|M| 



\ 



/ 



/ 



\ ,/\M\ J 
\ 



\ 



J 



a b 
6t c 



\" '/ 



l-\M\ 

/r^PT 



/ 



/ 



whose positivity follows immediately from the 2-positivity (in fact, complete positiv- 
ity) of the trace operator. 

For our purposes, we will choose p = pb(£) and a = PBj(e.) ® PsEii) &^d take M 
to be P'b, the non-trivial part of Pb = P'b ® 'lv\B{e), noting that ||P^|| < 1, since 
\\Pb\\ < 1. Letting pi = Tr{\P'j^\pB(i)) and p2 = Tr(|P^|pB^(£) ® pBsie)), we get from 



(3.3.22) and the definition of $: 



S {pB(i) WPBiii) ® Pbew) > -H2{pi) - pi In /i2 - (1 - Pi) ln(l - P2) > -pi In p2 - In 2, 

(3.3.24) 
where we used again the bound H2{c) < ln2 for the binary entropy and the obser- 
vation that (1 — yUi)ln(l — P2) < 0. We will now show that /ii > (1 — e(^))^ and 
P2 < (2P + 4e(£))/(l-e(£))2. We start with the lower bound for pi = Tr(|P^|pBw) = 
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Tr(|Pi,|Po) > Tr(|PB|2Po) = \\PoPb\\1 and hence, 

^/^^l > WPoPsh > WPoPEPyPywh > \\Poh " WPoiPEPyPyw - Po)h > 1 " ^W, 

(3.3.25) 
which translates, for i > ii such that e{i) < 1, into the lower bound /ii > (1 — e(^))^- 
For /i2, we note that Tr (|P^|pb^(^) ® PBEd)) = TrdPelpy ® Pv\y)- Moreover, given 
the polar decomposition Pb = U\Pb\, we have that 

Tt{\Pb\PyPv\yPy ® Pv\y) = 1:t{U^ PbPyPv\yPy ® Pv\y) 

< \\U^\\ I 1i{PbPyPv\yPy ® Pv\y)\ 

< Tr(Popy ® pv\y) + \\PbPyPv\y - Poll 

< P + e{i). 

Now, observe that setting /is = Tr(PyPv\ypy Pv\y) we have: 

P3 > I Tr(PyPy\yPo)p > | Tr (P^Py Py\y Pq) |' > (1 - t{^)f (3.3.26) 

from which we get: 

Tr [\Pb\PyPv\yPy ® Pv\y] = Tr [(|Pb| - P2) {PyPv\y - Ps) Py ® Pv\y] + P2P3- 

(3.3.27) 
Hence, by the triangle inequality, 

P2P3 < Tr [\Pb\PyPv\yPy ® Pv\y] + |Tr [(|Pij| - ps) {PyPv\y - Ps) Py ® Pv\y] \ ■ 

(3.3.28) 
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Moreover, by Cauchy-Schwarz, we have the following upper bound 



|Tr [(|Pb| - /is) {PyPv\y - /^s) Py ® Pv\y] 



< JTt [{\Pb\- 1^2) pY ® Pv\y] ■ A/ Tr {PyPv\y - Ps) Py ® Pv\y 



'Tt [\Pb\^Py ® Pv\y] - pI ■ \J Pz - p\ 
< ^Jii2 - pI ■ yP3 - pI< VM^y/P^ 



(3.3.29) 



where we used ||Pb|| < 1 in the previous to last inequality and 1 — P3 < 2e(£) in the 
last inequality. Setting x = ^^2Pz, the above inequalities imply: 



X 



- ^fMt)x - (P + e(£)) < 0, 



(3.3.30) 



and solving the quadratic inequality implies that x^ < 2P + 4e(£), and hence, 



P2 < 



2P + 4e(£) 

(i-^W)^' 



(3.3.31) 



Plugging things back in (3.3.24), yields (3.3.19). We are now ready to make use of 



bound (3.3.21). Setting 



Se = max|S'(pB^(^j(£)) : BY{m){^) C Py(„*/2)(m72) 



we get from (3.3.21): 



S21. < 2S'£ 



and one can easily check that for £2 < ^ < m* /2 we have: 



Sf< 



Se, Jg(7(4£2)) i2 



2^ 



e 



(3.3.32) 



(3.3.33) 
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Finally, since Sm'/2 > and S^^ < (i2'^YiN\dA\, the above bound implies that: 

m*<8£2 22«^°^l^^l =^ 5(p^) <32£2lnAr22«i°^l^^l|5A|, (3.3.34) 

which ultimately yields: 

S{pa) <32e(ln(12C(7,rf,J)) + 21n|M|)lnA^22«i"^l^^l|M|. (3.3.35) 

Obviously, the above bound gives a constant upper bound in one dimension. Never- 
theless, the exponential factor prohibits any direct extension of the area law to higher 
dimensions. D 

3.3.1 Proof of Entropy Bound 

In the proof of Case I in the above Theorem the following lemma was used. 

Lemma 3.3.2 (Entropy bound). Let p he a density matrix with spectral decom- 
position X]a'^('^)lV^('^))(V^('^)l ^''^'^ assume that there zs < c < 1, such that for an 
increasing sequence s„, n > 1 with s„+i/s„ < i?, Vn > 1 and si > 1, the following 
constraint holds on the eigenvalues cr{a): 



J2 ^(") < c", Vn > 1. (3.3.36) 

a>s„+l 

Then, the entropy of p is bounded by: 



S{p) <\ns^ + -^\nR+ -^H^il - c), 
1 — c 1 — c 



where H2{1 — c) = —cine — (1 — c) ln(l — c) is the binary entropy. 
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Proof. Since the Shannon entropy is Schur-concave, it is maximized when the input 
probabihty distribution {cr(a)} is "minimized per entry" (i.e. majorized by all other 



probability vectors satisfying the constraint ( 3.3.36[ )). Note that s„+i < Si-R", n > 1 
and hence s„+i — s„ < Si-R". 

For a < Si we have that — ^^ a{a) In o"(a;) is minimized when all a{a) are equal 
and their sum is minimized according to (3.3.36[): 



gi -, _ 
2_^o-o{a) = 1 — c =^ a{a) = . 

Similarly, for a G [s„ + 1, s„+i], ra > 1, we see that the entropy is minimized when 

J2 cx{a) = c" - c"+i =^ a{a) = ^ _^ . 

a=Sn + l 

Gathering terms, the entropy can be written as S{p) = J2'^=o^'n.{p), where S'o(p) = 
— ^^^^^^ cr(a;) lncr(a;) and Sn{p) = ~ J2'^a=sn+i '^ i^) ^^ '^ i^) ^ n > 1. We can bound 
each Snip) ^^ follows: 

Soip) < -(1 - c) ln(l - c) + (1 - c) In Si 



and for n > 1, 



Snip) < (1 - C)C^{- ln((l - C)cn + HSn+1 ' S„)) 

< c"(l - c)(lnsi - ln(l - c)) + nc"(l - c)(lni? - Inc) 
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Hence, we get the desired bound: 

oo 

S{p) < In Si - ln(l - c) + (1 - c)('^nc") (In i?- In c) 

n=l 

= lnsi + -^lnR+- H2{l-c). 

1 — c 1 — c 

and this completes the proof of this Lemma. D 

3.4 Approximation for the ground state of a gapped 
Hamiltonian 

In this section, we construct an approximation to the ground state of a gapped 
Hamiltonian. The setup is the same as in the previous section, but we restate the 
main assumptions here for the sake of completeness. Again, we consider a system of 
the following type: Let l^ be a finite subset of Z^. At each x & V, we have a finite- 
dimensional Hilbert space of dimension rix- Let Ax be the algebra of observables 
associated with X G V. We consider the dynamics r^ generated by the Hamiltonian 

xcv 

The main assumptions are as follows: The interactions ^{X) are uniformly bounded 
and finite range, and, for convenience, we will take pair interactions with range 1. 
So, for all X C V, $(X) = <I>(X)* G Ax, ||$x|| < J, for some constant J > 0, and 
$(X) = if diam(X) > 1. Here, for X cV, 

diam(X) = m.ax{d{x,y) \x,y E X}. 
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Moreover, the boundary of X, denoted by dX, is 

dA = {x E A \ there exists y E V \A, with d{x, y) < 1}, 

and we define the distance between sets X and Y to be 

d{X, Y) = mm{d{x, y) \x e X,y e Y}. 

For an observable A, we define t^{A) = q^^Ha j^^-ttHA ^ where H/^ = '^xca^(^)- 

In order to construct the approximation to the ground state given below, we will 
assume further that the Hamiltonian Hy has a unique, normalized ground state, 
which we will denote by |^o)) and a spectral gap 7 > to the first excited state. 

Theorem 3.4.1. There exists .^ > 0, such that for any sufficiently large i > c^d'^C,'^ 
there exist two orthogonal projections Pa G Aa, Pv\a £ Av\a o-nd an operator Pb G 
•AB{A;i) with II -Pb II < 1, such that 

\\PbPaPv\a- P4<C{l,d,J)\dA\''e-'/^ (3.4.1) 

where Pq is the projection onto the (unique) ground state, Cq is of order unity and 
C{'~f,d,J), C, are explicit in terms of J, 7, and the dimensionality ofV. 

The following lemma is an essential tool for the approximation theorem, so we 
state it and prove it here. 

Lemma 3.4.2 (Lieb-Robinson Bound). Let Hy be a Hamiltonian satisfying the above 
assumptions. Then, for a region A G V and two observables A G Ax and B G Ay, 
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with X C A and F fl A 7^ 0, we have the following bound: 

\\[rf'{A),B]\\ < 2p||||fi|||aX|e-^^(^'^)e"l*l, V|t| < e'^^^+^^ '^^^'^K (3.4.2) 

where v = 4{2d — 1)J is the Lieb-Robinson velocity of the system, and fi is a free 
parameter. 

We should note here that similar to the bound given in |16], the above bound is 
actually true for all times t, since for large |t| and fixed d{X, Y), fi becomes negative, 
making the above upper bound non-optimal, as one always has the naive bound 
||[r/^(y4), i?]|| < 2||74||||i?||. The above "constraint" in |t| is given only to clarify the 
regime in which the Lieb-Robinson bound can be more useful than the naive upper 
bound given above. 



Proof. Following closely the proof in |46] , we define the quantity 

Cb{Z; t) = sup 



[rnA),B]\ 



AeAz 



\A\ 



where B G Ay is fixed, and Z C A is considered arbitrary. For A G Ax, we 
introduce the operator function fit) = [tI^{t^^{A)),B]. Noting that r^tiA) G Ax 
and that t^^{[A,B]) = [Tf'{A),Tt^{B)] we may compute 



nt) = ^[[rf{Hs^x)),r,Hr''M))lB] 

= z [n^ (Hs^x)) , fit)] - ^ [rHr^M)), [^^ (^5(x)) ,B]], 

where we used the Jacobi identity [A, [B, C]] + [B, [C, A]] + [C, [A, B]] = for the 
second equality, and defined S'(X) to be the set of nearest-neighbor interactions 
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{Z G A\X : Z n dX 7^ 0}. Now, one may easily verify that f{t) is given by: 



e^-^W (^/(O) -^y^ e^^(-^) {[r^ir^AA)), [r^ {Hs^x)) ,B]]) e-^^^-^Us ) e^^-^W, 
where A{t) = J^ t^{Hs(x)) ds. Since A{t) is self-adjoint, we have that: 



ll/WII < 11/(0)11+/ \\[rt{r^s{A))Art{Hs^x)).B]]\\ds (3.4.3) 

< \\[A,B]\\+2\\A\\ J2 f\\[r^iHZ)),B]]\\ds (3.4.4) 



zesix)-^^ 



From (3.4.4) and the fact that ||r^^(y4)|| = ||v4||, we get the following recursive rela- 
tion: 

CB{X;t)<CB{X;0) + 2 ^ \mZ)\\ f CB{Z,s)ds (3.4.5) 

zes{x) -^° 

Moreover, it should be clear that Cb{Z, 0) < 2\\B\\6y{Z), where Sy{Z) = OiiYnZ = 



and 6y{Z) = 1, otherwise. Iterating (3.4.5) we get: 



CB(X;t)<2||5||V^^a„, (3.4.6) 



n=0 



where 

<^n= Y. E ■■■ E 5Y{Zn)Tiumz{)i 

Zi&3{X) Z2&SiZi) Z„eS(Zn_i) 

for n > 1. Since |5'(Z)| < 2{2d — 1) for Z with diam(Z) < 1, we have the following 
upper bound on a„: 

an < ({§^) [2{2d-l)Jr6y{S^-\X)), (3.4.7) 

where S^"\X) = S{S{- ■ ■ S{S{X)) ■■■)), composition n times. Now, note that the set 
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5'(")(^X) contains sites at most distance n away from X; that is d{S^"'\X),X) = n. 
Hence, d{S^''\X)X) > d{XX) - n, which imphes that a„ = 0, for n < d{X,Y)- 
Going back to (3.4.6) and noting that ( \^_2 ) < \dX\ we get: 



CB{X-t) < 2\\B\\\dX\ J2 



mt\y 



n\ 



< 2||E|||(9X| 

< nB\\\dx\ 



n=d(X,Y) 

(4(2rf- l)J|t|)'^(^'^) ^ (4(2rf- l)J|t|)' 



d{X,YV- 
d{X,YV. 



ra=0 



n! 



(3.4.8) 



where v = 4{2d — 1)J. Now, using the simple bound Inn! > nln(n/e) for n > 1, 
which one may easily prove by induction, we have that: 



d{x,Yy- - 



\t\ < e-(i+^) 



-(i+,)d{X,Y) 



(3.4.9) 



which combined with (3.4.8) and the definition of Cb{X; t) completes the proof. D 



We are now ready to prove Theorem 3.4.1 The proof consists of three main steps 



which translate into Propositions |3.4.4[ |3.4.6| and |3.4.8[ respectively. Here is a brief 
overview of the proof: 

The first step is to use the spectral gap in order to transform the original Hamil- 
tonian into one consisting of 3 parts, each of which has exponentially small energy in 
the ground state with respect to some length L In particular, the 3 pieces are chosen 
so that the middle piece, which is a boundary-like region of increasing thickness £, 
acts as an interaction bridge between the other two pieces of the Hamiltonian. 

Nevertheless, at this point, the transformation of the Hamiltonian will introduce 
long-range interactions within each of the 3 components, which becomes an obstacle 
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when one wishes to isolate the contribution of each piece to the ground state of the 
original Hamiltonian. In order to reclaim a strictly local range of interactions for each 



component, we use the Lieb-Robinson bound given in Lemma 3.4.2 to concentrate the 
interactions within each component in regions of larger, but strictly local, support. 

At this stage, we will have transformed the original Hamiltonian into a new one 
that is close in norm to the original and has, furthermore, the desirable property of 
consisting of components each of which has small energy in the ground state and 
whose support is localized within a region of our choice. It remains now to construct 
projectors onto the low energy states of each component and combine them to produce 
the desired approximation to the ground state. 

After this brief overview, we are now ready to prove the theorem. 

Proof. First we note that Hy can be partitioned as a sum 



Hy = Hj + Hb + He, 



(3.4.10) 



where 



and 



Hj= J2 "^W' 



XCV: 



Hb=Y. $(X), 



XCV: 

XCB 



he= Y. "^w- 



XCV: 



Moreover, we have that 



(3.4.11) 



(3.4.12) 



(3.4.13) 



Lemma 3.4.3. The following commutators have norms bounded by the number of 
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interactions at the boundary of the relevant regions: 

\\[Hv,Hi]\\ < 8d^J'\dI\, \\[Hv,He]\\ < 8d'j'\dE\, (3.4.14) 

and 

\\[Hv,Hb]\\ <8d^JWdI\ + \dE\). (3.4.15) 



We prove this lemma after the proof of Theorem |3.4.1 

Without loss of generality, we may add a constant to each of the above Hamil- 
tonians and thus assume that their ground state energy is zero. In other words, we 
really work with the Hamiltonians Hx = Hx — {^o-, Hx'^o), for X G {V, I, B, E}. 

In this case, each term has zero ground state expectation. However, it is not clear 
that these terms have a quantifiably small norm when applied to the ground state. 
To achieve small norms, we introduce non-local versions of these observables. Given 
any self-adjoint Hamiltonian H, any local observable O, and any a > 0, we may 
define the (non-local) observable 



a -"^ 



iO)a = d-J n{0)e-'''dt, (3.4.16) 

where here Tt{0) = e**^Oe~**'^. It is easy to see that for every a > 0, ||(0)a|| < ||0||. 
With respect to the Hamiltonian Hy and a > (a free parameter to be chosen 
later), we now consider the non-local operator 

{Hx)a = J^j^ n{Hx)e-'^''dt, (3.4.17) 

forXG {I,B,E}. 
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Note that these non-local observables still sum to the total Hamiltonian, i.e., 

{Hj)^ + {Hb)^ + {HE)a = Hy. (3.4.18) 

Also, they have zero ground state expectation, i.e., 

(V'O, (i^/)aV^0) = (V-O, (i^B)a^o) = (V'O, (i^E)aV'o) = 0. (3.4.19) 

Lastly, each of these observables, when applied to the ground state, has a small vector 
norm. 

Proposition 3.4.4. Under the assumptions above, we have that for any a > 0, 

||(i^x).^o|| < M^^^^^e-i^, (3.4.20) 

7 

forXe{I,B,E}. 



We will prove Proposition 3.4.4 after we finish the proof of this theorem. We want 
the above bound to have an exponential decay of the form e^^^^ for some positive ^' 
that will be determined later. Thus, we set 

a = ^. (3.4.21) 



This completes the first step in the proof of Theorem 3.4.1 The next step is to 
approximate each of {Hj)a, {HB)a, and {HE)a with observables M7(a),Ms(a) and 
Me^o). These observables will be constructed in such a way that each has support 
on a specific subset of V and a small vector norm when applied to the ground state. 
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For X G {I,B,E}, set 



a '°" 



Mxia) = d-J < (Hx) e^"* dt. (3.4.22) 

Here, for any local observable O, we have introduced the following evolutions 

al = Tf{0) = e^tHAQ^-itHA (3.4.23) 

af = rf^^'^^^(0) = e^*^s(^^2^)Oe-^*^s(^^2^) (3.4.24) 

af = r^\'^{^0) = e''"^\^Oe-''"^\^. (3.4.25) 

Observe that the local Hamiltonians defined above have been chosen such that 

supp(M/(a)) = A, supp(MB(a)) = B{A; 2i), and supp(M£;(a)) = V\A. 

We now state and prove an upper bound on the number of lattice points on the 
surface of a sphere in Z'^. 

Lemma 3.4.5. Let s{n,d) denote the number of points on a sphere of radius n in 
if- . Then, the following bound holds: 

s{n,d)<2V-\ (3.4.26) 

Proof. Spheres in Z*^ look like rf— dimensional diamonds. The sphere of radius n may 
be constructed by choosing any of the d hyperplanes going through an "equator" of 
a sphere of radius n — 1, and lifting the two "hemispheres" in a direction normal to 
the chosen hyperplane to insert an equator of radius n between them. Noting that 
an equator of radius n of a rf-dimensional sphere is just a (rf — l)-dimensional sphere 
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of radius n, we get the following recursive relation: 

s{n, d) = s{n-l, d)+s{n, d-l)+s{n-l, d-1), s{l, d) = 2d, s{n, 1) = 2, (3.4.27) 

where the term s{n — 1, c? — 1) comes from adding the missing equator to one of the 
two lifted hemispheres that make up the new polar caps. Note that the above bound 



is true for d = 1,2, with equality. Moreover, we have from (3.4.27) that 

n 

s{n,d) <s{n-l,d) + 2s{n,d-l) =^ s{n,d) <2'^s{k,d-l) + s{l,d)-2s{l,d-l), 

fc=i 

and noting that s(l, d) — 2s(l, d — 1) = —2{d — 2) < 0, for d > 2, we get 

n 

s{n, rf) < 2 ^ s{k, d-1) <2 ns{n, d - 1) < (2n)^"^ s{n, 1), d>2. 

k=l 



which proves (3.4.5). D 



We will now use the finite (Lieb- Robinson) velocity v = A{2d~ 1)J of our dynam- 
ics, to show that our newly defined local Hamiltonian operators are close in norm to 
the original ones. 

Proposition 3.4.6. Set ^' = A + A (-j . For X G {/, E}, the following bound holds 

\\{Hx)a - Mx{a)\\< 

e-^/«' { I -^==f/'\dX\+3''-'H''d''-'i''-'\dA\] ^^ S^ (3.4.28) 
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and similarly, 



iHB)a - MB{a)\\ < 



2 _L ^,2 



^V 7 + "^ 



f/WdB{A;i)\) + 3'^+H'^d'^-H'^-^\dB{A-2i)\ 



Ad'j- 



2 72 



The proof of Proposition 3.4.6 appears after the proof of Proposition 3.4.4 below 



Now, for any X G {/, 5, -E}, the bound 



\Mx{a)M\ < UHx)aM + UMx{a)-{Hx)a)^o\ 



(3.4.29) 



is trivial Thus, given Proposition 3.4.4 and 3.4.6 when apphed to the ground state, 
each of these local observables will have small vector norm along the perscribed 
parametrization. 



As an immediate consequence of Proposition 3.4.6 and Lemma 3.4.5 , we have the 
following corollary. 

Corollary 3.4.7. Set D = max{i,2(d - 1)} and 



Ci{^,d,J) 



2d+i 



tta/t^T^ 7 



_|_ ^ _|_ od+2 A2d rd-1 \ (2ttJ) 



Then, the bounds 



\Hv - {Mi{a) + Msia) + ME{a))\\ <2Ci{^,d,J)\dA\i^e-e (3.4.30) 



and 



max ||Mx(a)^o|| < Cih,d, J)\dA\re-^ 

X&{I,B,E} 



(3.4.31) 



are valid. 
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Proof. From the definitions of the relevant sets, it is clear that since for each point 



in dA there are s{i, d) sites at a distance i > 1, using Lemma 3.4.5 we have 



max{2\dA\,\dI\ + \dE\,\dB{A;i)\} < 2^r-'\dA\. 



(3.4.32) 



and similarly, 



\dB{A;2£)\ < 2 



2d-lnd-l 



\dA\. 



(3.4.33) 



Now, the first bound follows from (3.4.18) and Proposition 3.4.6, Similarly, the second 



bound is obtained by combining (3.4.29) with Propositions 3.4.4 and 3.4.6 Note that 



the constant Ci(7,(i, J) was chosen optimally for the second bound, but one could 



eliminate the — term for the first bound. 



D 



The next step in the proof of this theorem is to define an approximate ground 
state projector. Consider 



^iHvt^-at^ rft. 



(3.4.34) 



Observe that for any functions / and g, 



if, {Pa - Po)9) = ^1^ e'"*' f e^''*^(/' Pi^9)dt - if, Pog) 



oo n^ poo 
a 



e*^*e-"* dtd{f,PEg) 



e 4. d{f,PEg), 



(3.4.35) 



which then readily yields that 



Pa — Pq 



< e 4a = e f 



(3.4.36) 
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We now define a first approximation of this ground state projector. Set 



P. 



J{Mi+MB+ME)t^-at^ ^^ 



(3.4.37) 



where we have dropped the dependence of Mj, Mb, and Me on a. Clearly, 



Pa — Pq 



< 



P - P 

^ a ^ a 



Pa — Pq 



(3.4.38) 



The final term above we have bounded in (3.4.36). To see a bound on the first term, 
we introduce the function 



Ft{X) = e 



i\{Mi+MB+ME)t ^i(l-\)Hvt 



(3.4.39) 



One easily calculates that 



F[{\) = -^te'^^^''+^'^+^'^^'{Hv 



- {Mi + MB + ME)}e 



i{l-X)Hvt 



(3.4.40) 



Thus clearly. 



P - P 

^ a ^ a 



< 



< 



— oo 
oo 



< 



2Ci(7,rf,J) 



F,{1) - FMW e-"-'" dt 

Hn - {Mi + Mb + Me) 
Hn - {Mi + Mb + Me)\\ 



|t|e-°* dt 



2 I .,2 



^V7 +'^ 



(3.4.41) 
(3.4.42) 
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along the given parametrization for a and i^'. Thus, (3.4.38) now yields 



Pa — Pq 



< 



^3 ^"^"^^ \dA\ £^+V2 + i| ^-m'_ 



(3.4.43) 



We now define the spectral projections Pa and Pv\a which appear in the statement 
of Theorem 13.4. 1[ 
Denote by 

For X G {/, E}, define Pa and Pv\a to be the spectral projection, corresponding to 
the self- adjoint matrices Mj and Me, respectively, onto those eigenvalues less than 



c. In light of Corollary 3.4.7, the following bounds: 



:1-Pa)^o\\ < -||M,^o|| < e 



-£/2€' 



(3.4.44) 



easily follow, and the same bound holds for 1 — Pv\a- To see that we can always find 
a non-zero Pa (and similarly Py\^), we simply observe that the above inequality for 
Pa can be written as 

(^oIPaI^o) > l-e-'/2«' (3.4.45) 

which implies that Pa is non-zero. Moreover, it is obvious that both Pa and Pv\a 
satisfy the required support assumption. 

These spectral projections may be inserted into our previous estimates as 



PoPaPvxa — Pq 



< 



(P^ - Po) PaPv\a\ + ||Po (1 - PaPv\a)\\ , (3.4.46) 



here we have dropped the dependence on A and L The first term above is estimated 
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as in (3.4.43) above. Since 



;i - PaPv\a) = ^ {(1 - Pa){1 + Pv\a) + (1 - Pv\a){1 + Pa)} , (3.4.47) 



it is clear that 



|Po(l - PaPv\a)\\ < \\Po{1-Pa)\\ + \\Po{l - Pv\a)\\ < ^e-'/^^' , (3.4.48) 



using (3.4.44). Therefore, we now have that 



PoPaPvxa — Pq 



< (l^l(^lM= \dA\ £^+1/2 + 3^ e-'/'^'. (3.4.49) 



We are now left with the task of finding a local operator Pb G Ab{A;3£) that 
approximates Pa- In order to do so let us write 



PqPaPvxa 



^i(Mi+MB+ME)t ^-i{Mi+ME)t ^i{Mi+ME)t p p g-"*^ df 



(3.4.50) 



Since the supports of Mj and Me are disjoint, it is clear that 



^^iM,+ME)tp^p^^^ _ p^p^^^ ^ ^ ^^^M,tp^ _ p^) (e^M.^p^^^ + p^^^) 



+ - {e'^'-'Pv\A - Pv\a) [e'''''PA + Pa) 



Moreover, we have that 



|(/,(e^^^^*-l)P^( 



[\e^^'-l)d{f,P^PAC 
Jo 



< c\t\\\f\\\\g 



(3.4.51) 
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and similarly for |(/, (e**^^* — l)Pv'\Afl')|- If we now define the following operator, 

then we have just demonstrated that 

\\PBia)PAPv\A-Po\\ 

< \\PB{a)PAPv\A- PaPAPv\A\\ + \\PaPAPv\A - Po\\ 

< f^^hhddl \dA\ £^+1/2 + 3^ e-^/^«' 

< C2{i,d, J)\dA\ e^'^e-''^^\ (3.4.53) 

where we have used the simple bound D < 2d — 3/2 and set 

C2(7, d, J) = + 3 

in the last line. 

Given our arguments above, we need only show that Pb{(^) may be approximated 
by a local operator Pb- 

Proposition 3.4.8. There exists a local operator Pb with support in B{A;3£) and 
\\Pb\\ < 1; for which 

2 /93a!, 



\PB-PB{a)\\ < Cs{^,d,J)\dA\'e'''e-r. (3.4.54) 
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where 



C3{l,d,J) 



4'^J 



2v 



Vi^ 



+ 4*^ + 



2v 



7^ + V 



;Ci{j,d,J) 



and Ci(7,(i, J) is given in Corollary 3.^.1 • 



The proof of Proposition 3.4.8|can be found at the end of this section. It is clear 



that once we have proven Proposition 3.4.8, the estimate (3.4.1) readily follows from 
the triangle inequality 



PbPaPv\a - Poll < \\PB{a)PAPv\A - Poll + \\Pb - Psia 



< iC2ii,d,J) + Csi^,d,J))\dA\'t''e-^ (3.4.55) 



after setting i' = 3i, since Pb G B{A]i') in the statement of theorem, and setting 



^ = 25/4^ = 25(1- 



'v\2\ 



since the final decay above is of the form e ^s' and we wish 



to cover the factor i^'^ by choosing a slower exponential decay for i' large enough. 



This completes the proof of Theorem 3.4.1 



D 



Now we prove the various Propositions used in the proof of the Theorem. 



Proof of Lemma 3.4-3 Define 



Ci = [Hv,Hj], Cb = [Hv,Hb], and Ce = [Hv,He]. 



(3.4.56) 



From the definitions of the sets /, B, and E, it is easy to see that Cj = [Hb,Hi], 
Ce = [Hb,He], and Cb = [Hi,Hb] + [He,Hb]. Clearly, for Z G {I,E}, we have 
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that 



[Hz,Hb] = Yl M^)'^b] 



XCV: 



d 



zedZ i=l 

d 



X X Xz^{z + a) X [$({z, z + e^}), <l>({z + Ci, u;})] 






■weB: 

d{z+ei,w)=l 



+ 5Z 5Z^^^^^~^*) 5Z [^{{z,z-ei}),^{{z-ei,w])\ 



zedZ i=l 



(2-ei,to)=l 



from which we get: 



[Hz,Hb]\\ <2{\dZ\d{2d-l){2J^)) 



(3.4.57) 



The bound 



is now clear. 



\[Hz,Hb]\\ <8d'J'\dZ 



Here is the proof of Proposition 3.4.4 



n 
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Proof of Proposition 3.4-4' We begin by observing that 



\\{Hx)a'^0\\ < -\\Hv{Hx)a'^0\\ 

7 

= -\\[Hv,{Hx)am\ 

7 

= -||(Cx)a^o||. (3.4.58) 

7 

For the first inequahty above, we used that {Hx)a'^o projects off the ground state, 
since (^ol(^x)al^o) =0. 

The value of ||(Cx)a^I/o|| can be estimated using the spectral theorem. For any 
vector /, one has that 



" / at- 



a '°° 



{f,iCx)a^o) = \l-l e-'^' {f,niCx)^o)dt 

e--t\f^e'tH,Cx^o)dt 
= J^l^ e-*' 1^ e^'^'dif, PECx^o)dt 
= j^ J^j^ e-"*'e^*^ dt d{f, PeCx^o) 
= / e-^d{f,PECx^o). (3.4.59) 

J-y 

In the third equality above we have introduced the notation Pe for the spectral 
projection corresponding to the self- adjoint operator Hy, and we also used that 
(\l'o, Cx'^o) = 0. The last equality is a basic result concerning Fourier transforms of 
gaussians (and re-scaling), (see e.g. (5.59) in |17]). We conclude then that 

\{f,{Cx)ai^o)\ < e-- 11/11 llCxIl, (3.4.60) 
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and with / = (Cx)a^o we find that 



\\{Cx)aM < e-^||C. 



X 



(3.4.61) 



Putting everything together, we have shown that 



1 ~^ 

7 7 



\[Hv,Hx]\\ _'A 



7 



(3.4.62) 



D 



Here is the proof of Proposition 3.4.6 



Proof of Proposition 3.4-6 ■ Since our method of bounding this difference is similar 



for X G {/,i?,i5}, we will provide the details only in the case of X = J. We begin 
by observing that 

\\{Hi)^ - Mj{a)\\ < J^j^ \\n{Hj) - rl{Hj)\\ e-"*'rft. (3.4.63) 

To bound the integral above, we introduce a parameter T > 0. For |t| > T, we use 
that 



"1*1 



\r,(H,)-Tf{Hj)\\ < I j-X^s{H,) - tHH,)) 

< 2\\[Hy,Hj]\\\t\. 



ds 



(3.4.64) 



From this, it follows readily that 



WiHi 



r^HM ^ 



'"^^^dt < ^ll[-^^'^^]IL-aT2 



\t\>T 



i/avr 



(3.4.65) 
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Again, we want the exponential decay to be of the form e~^/^ , which imphes that 

2P 
T = — . (3.4.66) 

For |t| < T, the estimate below is an immediate consequence of Proposition 1 in 



M 



\\n{Hi) - Tf{Hi)\\ < [ \\[Hv-HA,rf{Hi)]\\ ds. 

Jo 

The above commutator may be written as 



(3.4.67) 



[Hy-HA,TfiHj)] = 5^5^XA=(a: + e,)[$({a;,x + e,}),r,^(iJ,) 

xi^dA i=l 

d 



x£dA i=l 



and for each x G dA, 



[<!>{{x,x±e,}),Tf{Hj)]= Y, [H{x,x±e,}),Tf {<!>{¥))]. (3.4.69) 



YCV: 



Each of these commutators we estimate using the Lieb- Robinson bound from Lemma 
remembering that it is useful for times |s| < e"'-^"'"'^-' £/v, since from the defini- 



3.4.2 



tion of /, we have d{Y, dA) > £ for F fl / 7^ 0. We may easily satisfy this requirement 



by choosing T < e ^^"^^•' i/v, which combined with (3.4.66) implies that 



Of 

T=^<e~^^+^U/v (3.4.70) 
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and, hence 



l/i' <e-^'+^^l/{2v) 



(3.4.71) 



Remembering that x G dA, Lemma 3.4.2 and (3.4.69) imply the following bound: 



d{y,y')=l 

< {2d)2\d{y,y'}\m{x,x±e,})\\m{y,y'}W\^\ 

j/e/ 

< 8cij2e-l^le-^^^e"''('^(^'^)~(^+^)), (3.4.72) 

ye/ 

For the last inequality, note that the number of sites y E I that are a distance 
(i(x, |/)=n + £+l, n>0 from x is bounded above by s{n + £ + 1, c?), the number 
of lattice points on a sphere of radius n + £ + 1 in 7/. We now use Lemma 3.4.5 to 
bound the sum ^^^^e-'^('^(^'S^)-(^+^» in (|3.4.72|): 



y&I 



-im 



n>0 

< 2"'£'^-^ ^(n + 2)'^-^e-^" 



n.>0 



< 2'^i'^-^e^'' 



d 



d-i 



9(-/i) 



d- 



vn>0 



.-M" 



where we used the bound i + n + 1 < i{n + 2), for i* > 1, ri > 0, for the second 
inequality. But, X]n>o^~'^" = (1 ~ e^'*)^^ and one can show inductively that 



^(1 - e-'^)-! <{d- 1)! (1 - e-^y^ < ^ ,'' e"'^ 



d{-lA 



/^ 
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This demonstrates that 



[mx, X ± e^}), rf (Hj)] \\ < 8d\ e"' (—\ J^ i<'~^ e'^'e''^'^ (3.4.73) 



and feeding back into (3.4.67) - (3.4.69), we get 



\niHj) - r,^(if/) II < 2d\dA\^^^ (—^ jH^-'e-^' (e^l*l - l) . (3.4.74) 



Letting C{d, jj) = 16 e^^ ( ^-^ ) d'^'^^ we may now estimate 
J^j_ \HHj) - rl{Hj)\\ e-^'"dt 



< ^!^i^ J' i'^-^ldAle-^' "" 



V 



,^^%-'''"dt 



-T 



< 






2 



-<'-^fdt 



C{d,fi) ^2 nd~l\^A\„~^ll^^^ 



V 



(3.4.75) 



At this point, since we are interested in exponential decay of the form e ^'^ , (3.4.21) 
gives the following constraint on ^': 



1 



/i 



l+(^ 



Combined with (3.4.71), the above constraint implies the following bound on /x: 



-'-^^' 



which is satisfied naively for /x = 1/4, since the r.h.s is bounded below by 1. Hence, 
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we finally have that 



e' = 4 1 + 



7 



(3.4.76) 



Adding our results for |t| > T and |t| < T, it is clear that along the given parametriza- 



tion (3.4.66) 



Hi)^ - MAa)\\ < 



2\\[Hv,Hi]\\ 3{12d) 



d+l 



'an 



+ 



rt-'\dA\\e 



^1^ 



(3.4.77) 



where we used the bound e'^/yU < 6, for /i = 1/4 to bound C{d,fi) in ( |3.4.75 ). The 

n 



bound claimed in (3.4.28) now follows from (3.4.21) and Lemma 3.4.3 



Here is the proof of Proposition 3.4.8 



Proof of Proposition 3.4-8. Our goal is to show that the operator 



Psia) 



^i(Mi+MB+ME)t ^-i{Mi+ME)t ^-at^ ^^ 



(3.4.78) 



is well approximated by a local observable; one with support in B{A; 3i). To do so, 
we will take the normalized partial trace over the complimentary Hilbert space, i.e., 
the one associated with B{A; 3iY C V. 

For our estimates, it is convenient to calculate the partial trace as an integral over 
the group of unitaries [8], see also [16]. Given an arbitrary observable A G Ay and a 
set Y G V, define 

(3.4.79) 



{A)y= / U*AUfi{dU), 

where Uiy^) denotes the group of unitary operators over the Hilbert space Hy^ and 
fi is the associated, normalized Haar measure. It is easy to see that for any A G Ay, 
the quantity {A)y has been localized to Y in the sense that {A)y G Ay. 
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Let us define 

Pb = {PBia))B(A;3i). (3.4.80) 

Clearly Pb has the desired support and, moreover, ||Pb|| < 1 since ||PB(a)|| < 1. 

The difference between Pb and PB{(y) may be written in terms of a commutator, 
i.e. as 

Pb - PB{a) = I U*[PB{a),U] fiidU). (3.4.81) 

Ju(B{A;3e)'=) 

Thus, to show that the difference between Pb and Pb{c() is small (in norm), we 
need only estimate the commutator of Pb{oi) with an arbitrary unitary supported in 
B{A-UY. 

This is easy to calculate. Note that 



a '°° 



V '' J — oo 

To estimate the integrand, we define the function 
A short calculation demonstrates that 

f(t) = % \f,mi^Ms^M^)ty^^^-i{M,^M^)t^jj-^ _ ^3^4^g4) 



The form of the derivative appearing in (3.4.84) suggests that we define the evolution 



at 



(A) = e^(*^^+*^s+A/B)t^g-i(Af,+MB+A^B)t^ for any local observable A. (3.4.85) 
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With this in mind, we rewrite 

= «ai(MB)/(t)+2[at(MB),f/]e'(*^^+^«+*^^)*e-'(^^+^^)*. (3.4.86) 

Written as above, the function / can be bounded using norm-preservation. In par- 
ticular, let Vii) be the unitary evolution that satisfies the time-dependent differential 
equation 

i—V{t) = V{t)atiMB) with V{0) = 1. (3.4.87) 

CLL 

Explicitly, one has that 

V(t) = (.i(Mi+ME)t^-i{Mi+MB+ME)t _ (3.4.88) 

Considering now the product 

gis) = Vis)fis) (3.4.89) 

it is easy to see that 

g'{s) = V'{s)f{s) + V{s)f'{s) = iV{s) [a^iMs), U] Vis)*. (3.4.90) 



Thus 

Vit)fit) = git) - giO) = [ g'is)ds, (3.4.91) 

Jo 

and therefore, 

fit) = iVit)* [ Vis)[asiMB),U]Vis)*ds. (3.4.92) 

Jo 
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The bound 



1*1 
< / \\[asiMB),U]\\ ds. 





(3.4.93) 



readily follows. 

Since the Lieb-Robinson velocity associated to the dynamics at(-) grows with 



we estimate (|3.4.93|) by comparing back to the original dynamics. Consider the 

hsir) = a,(r,_,(MB)), (3.4.94) 



interpolating dynamics 



for < r < s and Tr{A) = e^^^^Ae *^v^. With s fixed, it is easy to calculate 



K{r) = tas-r {[{Mj + Mb + Me) - Hy, Tr{MB)]) . (3.4.95) 



We conclude then that 



MMb) - Ts{Mb) 



h'g{r) dr 

< 2\\Mb\\\\Hv-{Mj + Mb + Me)\\s 

< 2{\B{A;2i)\J) (2Ci{-f,d,J)\dA\i^e-^/^') s 

< 4'^+VCi(7,(^,J)|M|2£^+V/«'s, (3.4.96) 



where for the second to last inequality above we used the definition of Mb and 



Corollary 3.4.7 and for the last inequality we used Lemma 3.4.5 to get \B{A;2i)\ < 
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{2i)2'^{2iY-^\dA\. Hence, 



r\t\ r\t\ 

||/(t)|| < / \\[r,{MB),U]\\ ds + WiasiMs) - t,{Mb),U]\\ ds 

Jo Jo 

r\t\ 
< / \\[r,{MB),U]\\ ds 
Jo 

+ 4^+VCi(7,(i,J)|M|2£^+''e-^/«'t2. (3.4.97) 



Plugging the above bound into (3.4.82), we find that 



\[PBia),U]\\ < 



"^^ dt 



< 



a 



■oo 
oo 



ltl 



,{Mb),U]\\ ds e-'^' dt 



-oo \^0 

Ci{-f,d,J)\dA\''i'^-e 



since 



oo 1 

t^e-'''"dt 

— oo 



2a' 



(3.4.98) 



We estimate the double integral above separately for |t| > T and \t\ < T, taking T 



according to (3.4.66). The integral for |t| > T is easily bounded as follows: 



1*1 



r, 



{Mb),U]\\ ds e""* dt 



< 2\\Mb\ 



\t\>T \J0 

a 



Itle""* dt 



\t\>T 

r)2d+l J 



'aTT 

22d+lj 



2 I ,,2 



T^yl + ^ 



\dA\l^^ 



1/2^-^/?' 



(3.4.99) 



To bound the integral for |t| < T we use Lemma 3.4.2 Remembering that supp(f/) C 
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B{A; 3iy and supp(Mb) C B{A; 2i) we have for |t| < e-^^+^''> 



\ts{Mb),U\\ < 2||MB||||f/|||a5(A;2£)|e-'^VI^I 
< 420|aA|2£2rf-ie-/^VI^I 



(3.4.100) 



Note that as we have aheady shown in the proof of Proposition 3.4.6, since |t| < T, 



the above bound is useful for /i = 1/4, and C,' given in the statement of Proposition 



3.4.6 Hence, 



J^J^ (/'*' ll[^«(^^)'^]ll ^^) ^""*'^^ 



< 4^^^ - \dA\^ fd-i^-fraciA^ /« / e^We-'^'^dt 



f 



J 



1 2 /j2d~l 



e 4e4a 



< A'"" - \dA\' 

V 
V 



(3.4.101) 



Combining (3.4.99) and (3.4.101) with (3.4.98) we get: 



[PB{a),U]\\ < 



4'^J 



2v 



+ 4« 



2v 



„2 I ^,2 



< C3(7,rf,J)|9A|2£3V^/«\ 



Ci(7,rf,J) |M|2^^+^+ie-^/«' 



where we set 



C3{l,d,J) 



4V 



2v 



tta/t^T" 



+ 4^ + 



2v 



7^ + 1>2 



Ci{i,d,J) 
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and used the bound D < 2d — 1 in the last hne. The estimate claimed in (3.4.54) 



now easily follows from (3.4.81). 



n 
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Chapter 4 

Multiplicativity of the maximal 
2-norm for depolarizing 
Werner-Holevo Channels 



4.1 Introduction 

In this chapteiQ we explore the question of multiphcativity for the output 2-norm of a 
class of quantum channels related closely to the c?- dimensional Werner-Holevo (WH) 
channel Wdip) = ^^((Tr(p))]l(i — p'^), which is known [62j to give a counterexample 
to the multiplicativity of the maximal output p-norm for p > 4.79, when d = 3. For 
large dimensions, the WH channel acts like the completely depolarizing channel, since 
ll^rf(p) ~ ^Irflloo = ^/d- This is a property shared by the recent counterexamples 
to multiplicativity for p > 2 by Winter and for 1 < p < 2, by Hayden [Ml l27] . 
Nevertheless, it has been shown [HI H] that VVd(p) satisfies multiplicativity for 1 < 



^Reprinted with permission from Journal of Mathematical Physics, 48, 122102. Copyright 2007, 
American Institute of Physics. 
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p<2. 

4.2 The setup and some useful lemmas 

It is a natural extension of the previous research on constructing counterexamples to 
multiplicativity, to study the output p-norm of channels of the form 

WA,d(p) = Ap+(l-A)W,(p), 

and ask if those channels satisfy multiplicativity for p- norms with p = 2. We adopt 
the name "depolarized" Werner-Holevo channels due to the close connection of VVa,^ 
with depolarizing channels in higher dimensions. We focus our attention to the study 
of the output 2-norm for the tensor product channel Wa,^ ® VVa,^ acting on bipartite 
states in Mct,{'C) ® Mrf(C) and show that multiplicativity is satisfied for this norm for 
all dimensions d. 

A direct computation of the eigenvalues of Wa,^ ® y^\,d{\'>Pi2){i^i2\) turns out to 
be much harder for < A < 1, than for the boundary cases A = 0, 1. The reason is 
that the output consists of a combination of the input state and its transpose/partial 
transpose, which in general do not share a common eigenbasis. To work around this 
difficulty, we compute exphcitly the output 2-norm of Wa,^ ® y^ \,di\''P 12) {''Pi2\) and 
the maximal output 2-norm of VVa,^ and study the difference 

VxMl2) = {\\WxA\lf - l|WA,d ® WA,d(|^12)(^12|)||2 (4.2.1) 

We show that "Z^A.d > for all input states and A G [0, 1], d > 2. We begin with the 
computation of ||WA,d||| in the following Lemma. 
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Lemma 4.2.1. The (squared) maximal output 2-norm ofWx^d is given by 

I|Wa,.L - — j^^ — 
Proof. It is easy to check that ||WA,d(|V')('^|)||2 is 

= TTiWx^mii'iY) 



< 



,2 , 2A(l-A)(l-|(^|^)p) (1-A)^ 
{d-2)\^ + l 



where \iIj) denotes the complex conjugate of {ip) in the standard basis. Taking {ip) = 

/!} , with |0), |1) two standard basis vectors, we see that equality can be achieved 

in the above expression and the result follows. D 

We now turn our attention to the more complicated output 2-norm of VVa,^ ® 

WA,d(|V'l2)(^12|). 

Lemma 4.2.2. The (squared) output 2-norm IIVVa,^® VVA,(i(|'?/'i2)(V^i2|)||2 ^■^ given by: 

+ 5^|(^12|«P 

- 2{S, + Rl){S, + {d-2)Ql){l-\\p4l) 

- S,\\Wx,41MpipI + P2P2), 

where Qx = ^, -Ra = A - Qx, Sx = 2\Qx, pi = Tr2|V^i2)(V^i2|, P2 = Tri|V^i2)(V^i2| 
and T denotes transposition. 
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Proof. It is easy to check that 

+ QxRx[pi(S)td + 'ii-d(S)P2\ 

+ Ql[td0td+Hh2){i^2\] 
- Y(l^l2)(^12ir 

where Ti,T2 denote partial transposition w.r.t. the 1*^*, 2'^'^ tensor factor, respectively. 
Taking the trace after squaring the above expression and noting that 

Tr |^i2)(^i2|(|^i2)(^i2ir'==Tr pfcp^, 

for k = 1,2 (which one can show using the Schmidt decomposition of \1p12)), we get 
the desired result. D 

The following observations will be very useful in the proof of the main theorem, 
so we state them here as lemmas. 

Lemma 4.2.3. Let ai < a2 < ■ ■ ■ < cTd be non-negative numbers that sum up to 1. 
Then, the following inequality holds: 

d 

2 



f^d > ^ CT, 



i=l 



Proof. The r.h.s. of the inequality can be thought of as the expected value of the 
random variable X given by Pr(X = cji) = ai. The upper bound then follows 
immediately. D 
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Lemma 4.2.4. Let ^ be a quantum channel and denote by $ the complex conjugate 
(w.r.t. the Kraus operators) channel. Then, the following inequality holds: 

11$ (g) $||2 < 11$ (g) $||2 

More importantly, the maximal output 2-norm on the r.h.s. is achieved on inputs 
with conjugate Schmidt bases: 

i 

where the conjugate is taken w.r.t. the standard basis. 

Proof. Let \il)i2) = Xli v^k«)l/«) be the Schmidt decomposition of the input. Com- 
puting 11$ ® $(|^i2)(^i2|)||2 and setting aij^k,i = Tr($(|ei)(ej-|)$(|efc)(Q|)), f3ij^k,i = 
Tr($(|/j)(/j|)$(|/fe)(/i|)) and o-jj-fc,/ = (Tiajakcri, we get: 

||$(g)$(|7/'i2)(V^12|)||2 = y^ \/<^i,j,k,l ai,j,k,l Pi,j,k,l 

i,j,k,l 

= / , \^^i,j,k,lT^e{aij^k,l Pi,j,k,l} 

i,j,k,l 

I |2 

^ V^ / \'^i,j,k,l\ 

- Z^ V(^i,j,k,l 2 

i,j,k,l 

I V^ / \f^i,j,k,l\ 

+ Z^ V(^i,j,k,l 2 

i,j,k,l 

— ^_^ V^^\jAi\^i,j,kA ) 
ij,k,l 

where we assumed the last inequality w.l.o.g. But, the last expression is equal to 
11$ ® ^{\'ipmax){'4'max\)\\2- Moreover, one can follow the above steps to show that for 
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every input |?/'i2) there is an input l^^ax) ^^^h that 

ii$®$(i^i2)(^i2i)ii^<ii$®$(i^:„...)(^Lj)ii^. 

This concludes the proof. D 

At this point, it is important to note that the above lemma and the fact that 

||$||2 = PII2 imply that ||<l> ® <l>||2 < ||$||i follows from ||$ ® $(|^max)(^ma:.|)||2 < 

ll'^lli- 

4.3 Proof of the Main Result 



In this section, we will show that the difference Vx^d defined in (4.2.1) is always non- 
negative, which is equivalent to multiplicativity of the output 2-norm for Wa,^. We 
state this as a theorem: 

Theorem 4.3.1. For the depolarized Werner- Holevo channel Wx.d, we have for A G 
[0,l],d>2: 

||WA,d®WA,d||2= IIWawII' 



Proof. From Lemma 4.2.2 we see that the condition ^A,d(|'^i2)) > is equivalent to 



Sl\{i^u\i^u)\' < 2(31 + P'x){l- Mil) 

+ SxWWxAlT^^iPiPJ + P2pI), 

where Pf = [Q| + {d - 2)Rl]Sx + {d - 2)Q^i?| > 0. Using Lemma 

II^A.dlli as (1 + y/d — l)Sx + (A — 7^)^) "^^ ^^^ ^^^^ it is sufficient to prove the 



4.2.1 



to write 
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following inequality 

\{^i2\M\' < 2(1 - llpii) + (1 + Vd^) TripipJ + P2pI) (4.3.1) 

(the boundary cases A = 0, 1 follow from 1 > ||pi||2)- 

Since both the identity and the WH channel remain unchanged under complex 
conjugation, their convex combination inherits that property. With this in mind, 



Lemma 4.2.4 implies that we only need to check if: 

|(^max|^/w)r < 2(1 - IIpiII^) + 2(1 + v^^^T) Tr(pipf), (4.3.2) 



since pi = pi for \i)rnax) = Z]iv^|ei)|ei). Moreover, ||pi||^ = ^i^^af, where some 



of the (Tj may be zero. Applying Lemma 4.2.3 (and borrowing its notation w.l.o.g.), 
it follows that ||pi||2 < cr^- It becomes clear now that in order to prove (4.3.1), it is 
sufficient to show: 

KV-maxlC^)!' < 2(1 - aa) + 2(1 + Vd^) Tr(pipf) (4.3.3) 



for (Td > 1/2, since \{4'max\'^max)\ < 1 and Tr(pipf) > 0. We now compute the 
following two quantities: 

\{^max\^max)\ = ^ y/o\0)j\{ei\e])\'^ (4.3.4) 

Tr(pipf) = $^v/^^Ke.|e7)l' (4-3.5) 



«j 



We will need to treat dimensions d < 4 and d > 5 separately. For c? < 4 we use 
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Cauchy-Schwarz to get the following estimate for 



(5Zy^^Kei|ei)P 



*j 



< (5^a,(T,|(e,|e-)|2)(^|(ei|e-)p 

< d'^aiaj\{ei\e])\'^ 



(4.3.6) 



where we have used Parseval's identity in the last inequality. Using (4.3.3) and (4.3.4) 



we see from estimate (4.3.6) that it is sufficient to show that d < 2(1 + \^d — 1) 



which is true for d < 4. We now turn our attention to the case d > 5. We will need 



a different estimate than the one given in (4.3.6), since we need to make use of the 



assumption that ad > 1/2 in order to lower the factor d in (4.3.6). We start by using 
Cauchy-Schwarz to get the following upper bound: 



where 



(5]^/^J^|(e.|e-)r)'<3(A2 + J22 + j|), 



«j 



h = ^y/^d^\{ed\ej)\'^ 

3 
i 
i^d,jj^d 



(4.3.7) 



A further application of Cauchy-Schwarz on /i = I2 and I3 will give us the desired 
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result. We start with an estimate for Ji. Noting that one of the summation indices 
is fixed to d, we get 

= ^(^d(Tj\{ed\e])\' 

j 
< Tr(pip[) 



|2 



Similarly, we have /| < Ti(pipJ). Since 1 + \^d — 1 > 3 for d > 5, we see from (4.3.3) 
and (4.3.7) that it remains to show 3/| < 2 (1 — aa)- We have 



, 2 

2 



^3 = [ Yl v^^\(^^\^j)\ 

s (E-.y 

= a-^d)' 



It remains to show that 3 (1 — cr^)^ < 2 (1 — cr^) <^ (1 — (Td)(3o"d — 1) > 0, which 
follows from our earlier observation that we only need to consider cr^; G [|, 1]. D 

4.4 Beyond entrywise positivity 

In this section, we show that the depolarized Werner-Holevo channels satisfy multi- 
plicativity despite having a more complex structure than quantum channels satisfying 
a condition sufficient for multiplicativity of the maximal output 2-norm. 
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Proposition 4.4.1. The depolarized Werner-Holevo channels VVa,^ with A G (0, 1) 
and d > 3 do not satisfy the following entrywise- positivity (EP) condition: There 
exists an orthonormal basis {\ei)}f^^ of C^ such that 

TTWx4\ei){ei\)Wx,d{\ej){ek\) >0, ^t,j,k,L 



Proof. One can check that Tt Wx^dH^i) {^i\)^x,d{\Gj) {^kl) is given by: 



A^ + 



■1-A\2i 



1-A\2i 



{-, SiiSkH -, - + (d-2){- 

\d-lJ '^ ' I d-1 ^ '\d-l 



2A(1 — A) . , — , ._, , 
^i,i^j,k T— I — {ei\ek){ei\ej) 



where \ek) denotes the complex conjugate of \ek), as before. Now, taking i = j,k j^ I 
in the above expression, we see that the EP condition imphes: 

{ei\e^){ei\ei) < 0, Wi,k ^ I. 

Summing over i in the above inequahty gives us 0, which imphes that: 



(ei|efc)(e/|ei) = 0, Vi, fc 7^ /. 



(4.4.1) 



Fixing /, we choose % = 7r(/) such that {ei\eT,(^i)) 7^ (we can always find such a 7r(/), 
since otherwise \ei) = 0; a contradiction to \ei) being an orthonormal basis vector). 



Condition (4.4.1) then implies that {e.,r{i)\ek) = 0, VA; 7^ /. Since the {|efc)} form 
an orthonormal basis, it follows that |e7r(/)) = |e/"),V/. We may now rewrite the EP 
condition as: 



A^ + 



1-A\2 



d-1 



^i,j^k,l + 



2A(1-A) 
d-1 



f(rf-2) 



1-A\2 



d-1 



Si,i5j,k > 



2A(1 - A) 
d-1 



Si,n{k)Sj,n{l) 
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Choosing i = 7i{k),j = 7r(/) and k ^ I, the above condition becomes: 
-2A(1-A) ,, ^/1-A\2]^ ^ 2A(1-A) 



d-1 



+ {d-2) f j 5^(fc),i5^(o,fc > 



d- 1 



The EP condition forces tt{1) = k^k ^ I (note that Tc{k) = I then follows from the 
definition of vr(A;),) which is impossible for d > 3. For d = 2, choosing |ei) = \^ 
satisfies the EP condition. D 



4.5 Discussion 

We have shown that for depolarized Werner- Holevo channels the maximum output 2- 
norm is multiplicative. For A G (0, 1) and d > 3, the depolarized Werner-Holevo maps 
do not satisfy the entry wise-positivity (EP) condition introduced by C King and 
M.B. Ruskai in BUBSl, which makes our result non-trivial. Moreover, a closer look at 



Lemma 4.2.4 in conjuction with random unitary channels may yield a counterexample 



for the multiplicativity of the maximal output 2-norm in the vein of 
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